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ABSTEACT 


This  thesis  is  the  first  study  in  a  bleached  pulp 
and  paper  mill  modelling  and  control  project  in  the 
Department  of  Chemical  Engineering,  University  of  Alberta. 

A  nonlinear  lumped  parameter  model  for  a  Kamyr  continuous 
kraft  pulping  digester  was  developed.  This  model  includes 
thirteen  state  variables  which  include:  the  wood  flew  rate, 
wood  component  concentrations  and  entrapped  and  free  liquor 
concentrations.  The  digester  is  split  up  into  a  number  of 
sections  each  of  which  are  an  ideal  stirred  tank.  These 
sections  are  further  divided  into  three  subsections  that 
include  diffusion  of  cooking  chemicals  and  reaction  products 
from  the  chip.  The  subsections  are  for  wood,  entrapped 
liquor  in  the  wood  and  the  free  liquor  flowing  in  the  chip 
void  volume.  The  nonlinear  set  of  ordinary  differential 
equations  are  solved  for  each  section  in  the  digester, 
including  the  digester  washing  sections. 

The  complete  digester  modelling  package  consists  of 
a  dynamic  and  steady  state  representation  of  a  continuous 
digester  and  one  of  a  batch  digester.  The  continuous 
digester  model  testing  was  done  in  parts.  Individual 
modules  of  the  program  were  tested  first,  followed  by 
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testing  of  the  overall  program  for  which  the  final  results 
were  known-  The  validity  of  the  kinetic  cooking  model  was 
compared  using  the  batch  digester  model  that  was  developed 
for  the  purpose.  The  dynamic  program  was  checked  by  running 
the  digester  from  starting  with  a  constant  initial  profile 
and  continuing  to  the  final  steady  state  and  comparing  this 
with  the  results  obtained  from  the  steady  state  model. 

Dynamic  open  loop  simulation  runs  for  the  Kamyr 
digester  were  conducted  with  the  ultimate  aim  of  evaluating 
various  industrial  control  practices.  The  three  main 
operating  control  variables  are  cooking  temperature,  cook 
zone  residence  time  and  cooking  chemical  concentration. 

These  variables  were  varied  to  evaluate  different  operating 
tradeoffs  and  the  effect  on  scalar  measures  that  are  used 
for  pulp  quality  characterization,  such  as  H-Pactor, 
residual  lignin  content.  Kappa  number.  The  H-Pactor,  for 
example,  that  incorporates  cooking  time  and  temperature  was 
found  to  be  inadequate  when  large  alkali  concentration 
changes  were  involved  and  the  residual  lignin  content  (Kappa 
number)  was  similar  although  significantly  different 
operating  conditions  were  used.  This  led  to  the  hypothesis 
that  "profile  control"  was  better  and  more  reliable  compared 
to  control  using  a  scalar  quantity.  The  implementation  of 
this  concept  should  be  practical  considering  the  current 
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availability  of  computer  hardware  and  technology. 

Open  loop  feedforward  type  of  control  was  studied 
for  a  chip  species  change  and  a  production  rate  change. 
Finally  typical  digester  closed  loop  control  schemes  are 
reviewed  and  presented. 
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CHAPTER  ONE 


INTRODUCTION 


A  pulp  and  paper  mill  modelling  and  control  study 
was  undertaken  to  provide  a  basis  for  further  work  in  the 
field  at  the  Department  of  Chemical  Engineering.  Pulp  and 
paper  production  is  an  important  segment  of  Canadian 
industry  and  it  was  felt  that  the  area  was  a  late  starter  in 
computer  control,  especially  when  compared  with  the 
petroleum  related  industries.  It  was  also  observed  that 
this  trend  was  bound  to  shift  in  the  future  as  faster 
minicomputers  and  better  sensors  became  available-  It  was 
also  obvious  that  in  the  meantime  a  lot  of  work  was  still 
required  to  develop  better  and  more  adaptable  models  and 
application  programs  of  various  pulping  processes  for 
computerized  control.  It  was  shown  in  earlier  studies  that 
the  economic  return  through  improved  regulatory  control 
justified  the  use  of  a  computer  (10,  43,  98).  This  project 
is  aimed  towards  a  study  of  some  of  the  underlying 
principles  and  problems  encountered  in  the  area  and  the 
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reasons  for  choosing  this  particular  topic  are  given  in  the 
following  paragraphs. 

The  digester  is  at  the  front  end  of  all  pulp  and 
paper  mill  unit  operations  (Figure  1.1)  and  therefore  is  a 
logical  place  to  start  a  pulp  and  paper  mill  project.  The 
pulp  from  the  digester  is  utilized  by  most  of  the  other 
sections  and  this  affects  the  performance  and  production 
capabilities  cf  all  downstream  units  right  down  to  the  paper 
machine,  which  makes  close  control  of  the  digester  all  the 
more  important.  It  is  a  reasonably  complicated  chemical 
engineering  process  that  involves  complex  physical  and 
chemical  relationships.  There  are  also  strong  interactions 
between  the  manipulated  variables,  and  the  digester  produces 
pure  time  delays  and  is  therefore  all  the  more  difficult  to 
control.  Operating  and  output  variables,  like  the  alkali 
and  wood  component  concentrations  are  unavailable  or 
unmeasurable  throughout  most  of  the  digester,  and  hence  the 
system  appears  very  much  like  a  *  black  box'.  A  systematic 
investigation  of  the  digester  should  lead  to  better  control 
by  the  computer. 

A  continuous  digester  was  selected  due  to  its 
advantages  over  a  batch  digester.  The  continuous  digester 
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is  at  a  disadvantage  when  production  flexibility  is 
required,  but  less  equipment  and  energy  is  reeded  for  the 
same  production  rate-  Nowadays  with  environmental  standards 
getting  stricter,  a  continuous  digester  performs  better  here 
also  because  of  the  relative  flow  rates  involved-  It  is 
estimated  that  a  500  air-dry-tons  (ADT) /  day  Kamyr 
continuous  digester  uses  124,000  fewer  barrels  of  fuel  oil 
per  year  when  compared  to  a  batch  digester/drum  washer 
system  (51).  The  pretreatment  of  chips  in  the  digester  feed 
system  gives  improved  liquor  penetration  and  the  high  heat 
washing  zone  reduces  the  load  on  downstream  washers.  This 
makes  it  possible  for  closer  control  and  these  days  the 
trend  is  definitely  towards  continuous  digesters  (51). 

Kraft  pulping  of  wood  was  chosen  as  it  is  tne  most 
dominant  in  chemical  pulping.  In  1974,  38%  of  the  total 
pulpwood  and  86%  cf  wood  residue  used  for  pulping  in  Canada 
was  for  Kraft  pulping,  this  is  equivalent  tc  56%  of  the 
total  production  (92)-  In  Sweden,  which  like  Canada  is  one 
of  the  major  pulp  and  paper  producers,  59%  of  all  pulp 
produced  in  1973  was  by  the  Kraft  process  (24  p134).  Kleppe 
(59)  has  a  table  for  1969  that  includes  wood  pulped  in  the 
United  States  and  it  shows  a  similar  breakdown- 
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1.1  Wood 

Wood  is  the  principal  raw  material  in  pulp  and 
paper  manufacture-  It  consists  of  fibres  joined  together  by 
a  cementing  material  and  numerous  minor  components  which 
vary  with  the  wood  species-  The  colourless  fibres  are 
carbohydrates  that  consist  of  relatively  long  chained 
cellulose  (5000-10000  monomer  units)  and  shorter  chained 
hemicellulose  (1000-2500  units).  The  coloured  binding 
material  is  lignin-  The  hemicelluloses ,  which  react  more 
rapidly  than  the  celluloses,  are  further  classified  as 
mannans  and  xylans,  whicn  are  short  forms  for  longer 
chemical  names.  The  extraneous  substances  include  terpenes, 
resins  and  phenols  which  are  soluble  in  organic  and 
inorganic  solvents  and  broadly  termed  as  extractives.  The 
main  objective  in  all  pulping  processes  is  to  separate  the 
wood  fibres  from  one  another  with  maximum  possible  yield  and 
minimum  less  of  fibre  strength-  These  processes  vary  from 
semi-mechanical  pulping  characterized  by  a  high  yield  of 
pulp  suitable  for  low  strength  products  such  as  fiberboard, 
to  comparatively  low  yield  chemical  pulping  for  fine  papers. 
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A  combination  of  mechanical  grinding  and  chemical  separation 
called  semi-chemical  are  between  these  two  extremes.  In 
chemical  pulping,  the  lignin  is  extracted  to  the  desired 
degree  by  using  either  an  alkaline  or  acidic  digesting 
liquor. 

1 . 2  Pu lp  and  Paper  Manufacture 

In  a  typical  bleached  pulp  and  paper  mill  ,  figure 
1. 1,  there  are  six  principal  distinct  steps; 

i)  Preparation  of  Raw  Materials, 

ii)  Pulping, 

iii)  Bleaching, 

iv)  Pulp  handling  and  or  Stock  preparation, 

v)  Paper  production, 

and  vi)  Recovery  of  chemicals. 

Logs  are  stored  in  the  millyard  before  deoarking 
and  chipping  in  the  woodroom.  The  chips  from  storage  are 
screened  to  remove  abnormal  sized  chips  and  undesirable 
material  and  then  moved  to  the  digester  where  they  are 
treated  with  cooking  chemicals  at  high  temperature  and 
pressure  to  produce  pulp.  Pulp  defibration  takes  place 
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Figure  1.1:  Kraft  Pulp  and  Paper  Manufacturing  Process 
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during  the  ’’blowing”  of  the  digester.  Pulp  from  the 
digester  is  then  sent  for  screening  operations  in  order  to 
remove  knots,  bark  fragments  and  uncooked  material  which  are 
recycled  back  tc  the  digester.  The  pulp  is  then  pumped  to 
the  washers  where  spent  chemicals  are  removed.  Pulp  washing 
is  performed  in  the  continuous  digester,  vacuum  filters  or 
in  diffuser  washers. 

Washed,  unbleached  pulp  is  used  for  kraft 
linerboard  and  bag  paper  whereas  pulps  that  require  a  high 
degree  of  brightness  are  subjected  to  bleaching  with  various 
chemicals.  The  bleaching  process  consists  of  a  number  of 
stages  of  reactors  and  vacuum  washers  where  the  pulp  is 
bleached  to  the  required  terminal  brightness.  Typical 
bleaching  sequences  are  CEHD,  CEDED,  and  CEHDED  where  C 
designates  chlorination,  E  is  for  caustic  extraction,  H  is 
for  sodium  hypochlorite  action  and  D,  Chlorine  Dioxide 
treatment . 

Before  paper  manufacture,  the  pulp  passes  through  a 
stock  preparation  stage.  This  involves  refining  and 
addition  of  chemicals.  The  refining  includes  beating  of  the 
pulp  whereby  the  wood  fibres  are  completely  separated  and 
chemicals  are  added  to  impart  desired  properties  to  the 
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paper  sheet.  The  pulp  slurry  is  then  passed  through  the 
headbox  onto  a  mesh  screen  to  remove  excess  water.  The 
resultant  fibre  mat  is  then  pressed  between  heated  rollers 
and  dried.  Finishing  operations  are  next  and  the  paper  is 
coated,  cut  to  size,  rolled,  weighed  and  moved  to  the 
warehouse  for  shipment. 

Except  for  some  small  amounts  of  chemicals  that  are 
not  recoverable  in  the  washing  stages,  the  major  part  of  the 
chemicals  used  in  the  digesters  are  recovered  and  returned 
to  the  process.  In  the  kraft  process  up  to  99%  of  the 
cooking  chemicals  can  he  recovered.  The  spent  chemicals  or 
the  black  liguor  is  evaporated  so  that  the  resultant  mass 
which  has  a  chemical  content  of  60-65%  can  be  burnt  in  the 
recovery  boiler.  The  smelt  from  the  boiler,  after  being 
dissolved  in  water  ,  called  the  green  liguor,  is 
recausticized  with  lime  to  make  fresh  cooking  liguor  or 
white  liguor. 


■ 


9 


1*3  Ka  myr  Digester 

The  Kamyr  digester.  Figure  1.2,  is  a  hollow, 
cylindrical  vertical  shell  where  wood  chips  and  the  cooking 
chemicals,  sodium  hydroxide  and  sodium  sulphide  react 
continuously  tc  produce  pulp.  The  chips  from  storage  are 
conveyed  into  the  chip  bin  and  then  into  a  hopper  which 
feeds  the  pockets  of  a  rotating  chipmeter  that  regulates 
flow  into  the  digester.  Chips  leaving  the  meter  are 
pretreated  in  the  steaming  vessel  before  falling  through  the 
chip  chute  into  the  high  pressure  feeder  that  sends  them  to 
the  top  of  the  digester.  In  the  initial  stages  the  chips 
absorb  the  cooking  liquor  without  much  reaction  nut  as  they 
travel  further  down  the  digester  the  temperature  is 
increased  ny  external  neating  in  two  steps.  This  is  done  by 
withdrawing  the  cooking  liquor  through  cylindrical  screens 
in  the  digester  periphery  and  passing  it  through  heat 
exchangers  before  returning  it  into  the  center  of  the 
digester,  slightly  above  the  screens.  At  a  temperature  of 
about  150  °C  the  rate  of  pulping  increases  considerably. 

The  cooking  time  at  this  temperature  depends  on  the  product 
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specif ications,  usually  the  yield  and  the  Kappa  Number,  and 
the  production  rate.  The  Kappa  Number  and  sometimes  the 
Permanganate  Number  are  used  as  different  measures  of 
residual  lignin  in  the  pulp.  The  cooked  chips  are  then 
quenched  by  injection  of  cool  liquor  to  stop  the  reaction 
and  the  spent  liquor  is  removed  at  the  extraction  screens 
below.  This  rate  of  liquor  removal  is  larger  than  the 
flowrate  in  the  cooking  stage  and  thus  results  in  a 
countercurrent  flow  of  dilute  liquor  in  the  wash  zone  lower 
down  in  the  digester.  The  spent  liquor  is  sent  to  the  flash 
tanks  and  evaporators  for  recovery.  The  countercurrent 
flowing  liquor  to  the  pulp  acts  as  a  washing  medium  and 
removes  reaction  products  and  unreacted  chemicals  from  it. 
Near  the  bottom  of  the  digester  the  wash  liquor  is  heated 
externally  to  improve  washing  efficiency.  The  pulp  is  then 
cooled  at  the  bottom  and  discharged  through  the  outlet 
device  and  blow  unit  to  the  blow  tank.  It  is  then  sent  for 
further  washing,  and  bleaching  if  so  required. 
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1.4  Ob  jecti ves  of  Study 

The  objectives  of  this  study  were  to  document  the 
state  of  the  art  and  current  operational  practices  for  a 
Kamyr  continuous  digester  and  to  derive  a  mathematical  model 
based  on  the  information  available.  This  model  was  to  be 
used  to  cneck  out  various  open  loop  schemes  and  to  work  out 
a  method  for  improving  digester  control  .  Part  of  the  work 
was  directed  towards  evaluating  schemes  such  as  "profile 
control",  and  multivariable  computer  control  techniques  that 
have  worked  well  in  other  applications  at  the  University  of 
Alberta.  It  was  hoped  that  simulation  studies  would  suggest 
an  improved  control  scheme  for  the  Kamyr  digester. 


CHAPTER  TWO 


DIGESTER  LITERATURE  REVIEW 


Alkaline  pulping  has  been  in  prominence  for  nearly 
a  hundred  years  now  (59,  81  p279)  and  during  this  period 
investigators  have  tried  to  explain  the  reaction  mechanism 
involved  in  wood  pulping  and  derive  applicable  kinetic 
models.  A  number  of  review  papers  have  been  published  on 
this  subject  in  recent  years.  The  Institute  of  Paper 
Chemistry  has  published  a  yearly  review  of  progress  in 
different  aspects  of  alkaline  pulping  from  1947  to  1968 
(46).  Other  papers  by  Wilder  and  Daleski  (101)  and  Kleppe 
(59)  provide  information  up  to  the  time  of  publication.  The 
textbooks  on  pulping  processes  by  Rydholm  (81)  and  "Pulp  and 
Paper  Manufacture"  (70)  have  an  extensive  list  of  references 
also. 
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2. 1  Review  cf  Continuous  Digesters 

The  first  Karoyr  digester  was  installed  in  Sweden  in 
1949  and  the  first  in  North  America  at  Hinton,  Alberta  eight 
years  later.  Smith  and  McCorry  (89)  gave  a  description  and 
operating  procedures  for  the  Hinton  digester.  Later  changes 
and  modifications  in  the  digester  cooking  system  are 
discussed  in  chronological  order  by  Jansson  (47) ,  starting 
from  the  hot  blow  of  pulp  to  the  cold  blow,  diffusion 
extraction  and  high  heat  diffusion  washing.  Knutsson  and 
Stockman  (60)  and  McKibbins  (74)  provide  evidence  to  support 
the  improvement  due  to  these  innovations.  ether  types  of 
continuous  digesters  are  the  steam  phase  digester,  a  two 
vessel  high  yield  system  with  a  separate  pre-impregnation 
vessel  and  the  saw  dust  digester  with  internal  washing  but 
the  emphasis  in  this  chapter  will  be  on  the  conventional 
Kamyr  hydraulic  digester. 

The  continuous  cooking  process  is  also  described  in 
various  texts  and  articles  (12,  49,  70,  85,  88,  93)  and  an 
operating  manual  on  this  digester  is  published  by  the  Scott 
Paper  Company  and  distributed  by  Kamyr  Inc  (85) .  Hubbe 
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2.2  digester  Modelling  Studies 

Continuous  digesters  have  been  modelled  with 
various  degrees  of  complexity  by  researchers  and  computer 
vendors.  Westinghouse  installed  one  of  the  first  digester 
computer  systems  at  Gulf  States  Paper  Corporation  (28)  but 
most  of  the  modelling  information  was  proprietary  and  not 
disclosed.  Smith  (90)  has  discussed  the  installation  and 
and  basic  formulation  of  this  model.  Johnson  (48)  derived  a 
linear  model  that  is  conceptually  the  same  as  that  at  Gulf 
States  though  it  has  less  detail.  He  developed  control 
equations  using  statistical  regression  methods  for  Kappa 
Number  with  cooking  temperature  and  chemical  concentrations 
as  the  dependent  variables.  His  analysis  was  based  on  batch 
cooking  data  augumented  with  random  noise.  Similar 


e 

. 


' 


. 


16 


statistical  work  has  also  been  done  by  Hoag  (40) ,  McKibbins 
et  al  (75) ,  Hatton  (36,  37)  and  Kerr  (52,  53)  and  empirical 
equations  developed  relating  the  important  digester  cooking 
variables.  Tucker  (98)  in  a  feasibility  study  for  process 
computer  control  mentions  a  related  model  to  the  earlier 
ones  for  quality  control  and  also  includes  a  measure  of  the 
time  delays  involved.  later  more  exhaustive  work  based  on 
chemical  engineering  principles  was  done  by  Stillman  (93) , 
Johnsson  (49)  and  Smith  and  Williams  (88)  and  will  re 
discussed  in  more  detail  in  the  next  section. 

The  digester  wash  zone  has  been  modelled  separately 
by  Williams  et  al  (103)  as  a  distributed  parameter  system 
and  steady  state  analytical  solutions  at  various  flow  rates 
are  given.  Separate  washing  studies  have  also  been  done  and 
a  dynamic  model  which  includes  cooking  in  the  wash  zone  has 
been  presented  by  Wetlessen  (100).  A  review  of  washing 
literature  and  a  survey  of  countercurrent  pulp  washing 
theory  and  models  has  been  presented  by  Toiriak  and  Lauzon 
(96,  97) .  They  have  an  extensive  list  of  references  and 
have  presented  the  three  basic  models  used  in  countercurrent 
washing  theory.  These  models  are  all  essentially  a 
variation  of  perfect  mixing  cells  in  series  with  the  washing 
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liquid  and  pulp  flowing  in  opposite  directions  to  each 
other. 


2.3  Survey  of  Past  Modelling  Work 

The  major  digester  models  developed  by  Stillman, 
Johnsson  and  Smith  are  discussed  in  detail  and  compared  in 
this  section.  Table  1  summarises  the  assumptions  and  main 
features  of  each  model. 

2.3.1  Stillman's  Model  (93) 


Stillman's  work  consists  of  two  parts.  A  batch 
digester  model  was  derived  first  to  test  the  cooking  model 
and  then  this  kinetic  model  was  used  to  derive  a  steady 
state  model  for  the  continuous  case.  The  batch  or  cooking 
model  consists  of  first  order  kinetic  rate  expressions  that 
were  derived  from  the  literature  existing  at  the  time. 

Seven  wood  components  were  considered,  for  the  broad 
classification  or  hardwood  and  softwood  species  namely, 
lignin,  cellulose,  mannan,  xylan,  acetyl  groups,  uronic  acid 
and  extractives.  The  reaction  rate  constants  were  derived 
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from  laboratory  cooking  data  and  these  constants  include  the 
effect  of  temperature  via  Arrhenius  type  relationships,  and 
the  effect  cf  varying  chemical  concentrations  in  the  cooking 
liquor  .  A  sodium  sulphide  hydrolysis  equation  is  used  to 
calculate  the  amount  of  alkali  instead  of  assuming  fifty 
percent  hydrolysis  to  the  hydroxide  and  a  factor  for 
evaluating  the  effective  sodium  sulphide-carbonydra te 
reaction  was  developed  after  fitting  cooking  data.  A  liquor 
density  relationship  is  also  derived  based  cn  the  expansion 
of  water  with  temperature-  The  batch  model  so  derived  was 
tested  against  punlished  results  of  hardwood  and  softwood 
cooks.  The  results  compared  fairly  well,  especially  in  tne 
final  stages  of  the  cook. 

JJhile  deriving  the  kinetic  constants,  Stillman 
mentions  that  although  substantial  literature  existed  on 
kraft  pulping  it  could  not  be  fully  utilised  because  of  the 
manner  in  which  it  was  collected.  The  main  problems  were  as 
a  result  of  lew  liquor  to  wood  ratios,  insufficient  cooking 
liquor  pre-impregnation  and  unavailability  of  liquor 
concentrations  during  the  cook  and  time  to  reach  the  maximum 
temperature . 

The  continuous  digester  model  consists  of  a  series 
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of  axial  sections.  Steady  state  material  and  energy 
.balances  are  conducted  on  each  section  where  wood  chips  and 
cooking  liquor  are  assumed  to  flow  with  the  same  velocity. 
Most  of  the  equations  derived  for  the  batch  model  are  used 
and  an  energy  balance  equation  was  the  only  extra  one 
derived.  The  cooking  equations  are  solved  over  the 
residence  time  of  the  pulp  in  each  of  the  axial  digester 
sections  instead  of  the  cooking  time  as  was  done  in  the 
batch  model.  The  model  testing  results  are  not  published  as 
customer  supplied  data  was  used  but  sample  calculations  for 
total  yield  based  on  Hubbe's  (43)  data  seem  to  agree  well. 

A  comparison  of  one  and  two  temperature  cooks  was 
made  and  a  higher  yield  claimed  for  the  same  lignin  content 
in  the  two  temperature  cook.  A  two  temperature  cook  is 
where  a  high  cooking  temperature  is  followed  by  a  lower  one 
for  the  remaining  portion  of  the  cook,  as  opposed  to  a 
higher  temperature  in  the  latter  portion  of  conventional 
continuous  cooks.  The  detrimental  effects  of  a  chip 
moisture  disturbance  on  pulp  production  rate  and  product 
quality  were  also  studied  by  simulation  from  this  model. 
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2.3*2  Johnsson*s  Model  (49) 

This  model  has  two  parts  to  it  also.  The  first  is 
a  detailed  distributed  parameter  model  for  all  sections  in 
the  digester  starting  from  the  steaming  vessel  to  the 
digester  outlet  valve.  This  model  includes  diffusion  of 
chemicals  and  reaction  products  into  the  wood  chip  and 
involves  a  measure  of  chip  thickness  distributions.  The 
second  part  is  where  the  former  dynamic  model  is  simplified 
to  a  steady  state  one.  In  this  static  model,  the  wash  zone 
Is  further  simplified  because  of  its  countercurrent  flow. 

Mass  and  energy  balances  are  derived  for  the  wood 
and  free  flowing  liguor.  The  wood  is  considered  to  have 
four  constituents,  lignin,  carbohydrates,  non  extractable 
solids  and  easily  extractable  solids  like  acetyls  and 
resins.  The  free  liquor  flows  in  the  chip  voids  and  the 
chemical  concentration  is  assumed  to  be  the  hydroxide 
effective  alkali  alone  since  the  effect  of  bisulphide 
concentration  is  neglected.  Effective  alkali  is  defined  as 
the  total  sodium  hydroxide  concentration,  assuming  that 
fifty  percent  of  the  sodium  sulphide  (Na2S)is  hydrolysed  to 
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sodium  hydroxide  (NaOH)  and  is  expressed  in  terms  of  sodium 
oxide  (Na20).  Reaction  kinetics  in  the  model  are  first 
order  in  lignin  and  the  carbohydrates  are  monitored  by  a 
lignin-yield  functional  relation ship,  that  is  assumed  known 
for  the  particular  wood  species. 

Partial  differential  eguations  for  the  model  are 
derived  with  lignin,  carbohydrates,  yield,  reaction  products 
and  effective  alkali  as  the  dependent  variables.  Tnese 
equations  are  similar  in  form  in  the  top  half  zones  of  the 
digester  with  the  cocurrent  flow,  as  the  only  difference  is 
with  respect  to  the  section's  temperature.  The  temperature 
profile  and  residence  time  relationship  are  assumed  known 
and  specified.  The  external  heater  streams  are  assumed  to 
be  perfectly  mixed  with  the  free  liquor  in  the  vicinity  of 
the  screens.  Eguations  for  the  wash  zone  are  derived  from 
those  for  the  cocking  zone  by  neglecting  the  pulping 
reactions  and  reversing  the  sign  of  the  appropriate  flow 
terms  to  account  for  the  countercurrent  free  flow  of  wash 
liquor.  This  simplifies  the  model  to  essentially  include 
only  the  diffusion  of  chemicals  and  reaction  products  out  of 
the  cooked  chips.  Johnsson  mentions  that  a  model  that 
includes  cooking  in  the  wash  zone  has  also  been  derived  but 
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was  not  included  in  the  published  work. 

The  dynamic  model  is  reduced  to  a  steady  state 
model  because  of  computational  convenience  and  for  testing 
and  design  purposes,  Johnsson  mentions  that  this  static 
model  is  fine  for  design  but  not  for  control  purposes.  The 
static  model  includes  the  basic  nonlinearities  of  the 
dynamic  model,  the  same  reaction  kinetics,  and  the  same 
dependence  cn  temperature  and  chip  thickness  distrit ution. 
Some  of  the  assumptions  made  to  obtain  a  solution  are  a 
constant  production  rate,  absence  of  disturbances  in  the 
input  variables  and  wood  parameters  and  a  constant  chip  size 
distribution  throughout  the  digester.  The  static  wash  zone 
problem  is  converted  from  a  boundary  value  problem  to  an 
initial  value  one  by  parameter  transformation. 

The  Johnsson  model  verification  and  testing  was 
done  using  data  from  a  high  yield  kraft  mill  digester 
operated  at  a  higher  and  lower  than  normal  production  rate. 
The  major  discrepancies  were  between  the  observed  and 
calculated  Kappa  numbers.  This  was  attributed  to  the 
presence  of  cooking  in  the  wash  zone,  errors  in  plant 
effective  alkali  measuremen ts,  diffusion  coefficients  and 
wood  density  values  used  and  approximations  in  the  model. 
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Several  studies  with  parameter  variation  and 
disturbances  in  the  input  variables  were  tried.  They 
include  variations  in  temperature,  lignin  and  acetyl 
content,  ligucr/wood  ratio,  alkalinity,  rate  coefficients, 
density  change  and  the  amount  of  void  volume  in  the  chip 
pack.  The  wash  studies  include  the  effect  of  washing 
variables,  such  as  dilution  factor  or  wash  flow  on  washing 
efficiencies. 

He  concludes  that  alkalinity,  ligucr/wood  ratio, 
temperature  and  cooking  time  affect  the  guality  and  cost  of 
production  of  the  pulp  and  that  chip  thickness  distributions 
are  helpful  in  studying  pulp  inhomogeneity  and  washing 
efficiency. 

2.3.3  Smith  and  Killian^s  Model  (88) 

This  is  a  dynamic  model  with  a  series  of  lumped, 
axial  sections  starting  from  the  digester  top  to  the  exit. 
Each  of  the  sections  in  the  model  is  assumed  to  include 
three  phases.  Therefore  a  section  is  further  divided  into 
three  subsections,  the  wood,  the  absorbed  or  entrapped 
liguor  in  the  wood,  and  the  cooking  liguor  between  tne  chip 
voids  that  is  free  to  move  about  and  is  termed  the  free 
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liquor.  The  model  considers  twelve  state  variables,  five 
for  the  wood  components,  three  each  for  the  component 
concentrations  in  the  free  and  the  entrapped  liquors,  and 
the  temperature  of  the  section.  The  five  wood  components 
considered  are  fast  and  slow  reacting  lignin,  and  the  three 
carbohydrate  fractions,  cellulose,  mannan  and  xylan.  The 
liquor  component  concentrations  included  are  hydroxide 
effective  alkali  and  bisulphide  ions,  from  the  hydrolysis  of 
sodium  sulphide  and  the  reaction  products.  The  temperature 
in  the  section  is  assumed  to  be  identical  in  the  three 
phases  (the  wood,  entrapped  and  free  liquors)  except  when  an 
externally  heated  stream  is  introduced.  The  external  heater 
stream  is  included  in  the  free  flowing  liguor  and  complete 
mixing  between  the  two  is  assumed.  The  dynamic  mass  and 
energy  balance  equations  derived  for  each  of  the  sections 
assume  that  it  is  an  ideal  stirred  tank  reactor  so  that  the 
exit  streams  have  the  same  conditions  as  that  in  the 
section.  They  claim  that  the  kinetic  model  that  was 
developed  is  based  on  observed  chemical  principles. 

Eeaction  rate  constants  are  developed  from  laboratory  and 
industrial  data  for  various  hardwood  and  softwood  species. 
Other  assumptions  made  in  deriving  the  model  are  that  the 
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entering  and  exit  flows  occur  in  their  respective  phases  and 
subsections  and  that  there  is  no  accumulation  of  wood  and 
entrapped  liquor,  this  gives  constant  flow  rates  to  the  wood 
and  entrapped  liquor  and  consequently  the  same  velocity  and 
residence  time.  The  residence  time  data  used  was  taken  from 
Johnsson  (49),  who  has  presented  a  curve  for  chip  volume  fed 
versus  digester  elevation.  The  free  liquor  rates  are 
constant  also,  except  where  there  are  external  flows  and  in 
that  case  they  are  taken  as  the  sum  total  of  both  flows. 

The  mass  transfer  rate  between  the  entrapped  and  free 
liquors  is  proportional  to  the  concentration  difference 
between  the  two  phases  and  the  rate  coefficients  are  assumed 
to  be  identical  for  all  the  components.  The  mass  transfer 
coefficients  have  a  temperature  dependence  and  piecewise 
curves  were  fitted  to  data  derived  from  the  Chilton  Colburn 
analogy. 

The  wash  zone  model  is  simplified  by  neglecting  the 
cooking  reaction  and  considering  only  the  dissolved  solid 
concentration  and  the  temperature.  This  reduces  the  number 
of  state  variables  in  the  model  to  three.  Dissolved  solids 
or  reaction  product  concentrations  in  the  entrapped  and  free 
liquors  and  the  section  temperatures  were  calculated  for 
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different  wash  flews  and  the  results  plotted  as  washing 
efriciency  and  wash  flow  rate.  Clean  wash  water  studies 
were  done  also  and  a  comparison  made  with  those  using  dilute 
liquor.  They  conclude  that  clean  water  washes  are  more 
economical  than  the  conventional  dilute  liquor  washes,  since 
tetter  washing  efficiencies  for  the  same  flew  are  obtained 
and  thus  save  energy  and  evaporation  costs. 

Testing  of  the  Smith  model  was  dene  with  actual 
industrial  data  and  the  results  agreed  fairly  well.  The 
kinetic  model  results  had  deviations  from  actual  laboratory 
data  in  some  cases,  hut  this  is  termed  to  be  the  best 
possible  fit.  Dynamic  simulations  were  done  by  varying  the 
chip  moisture,  alkali  content,  the  lower  heater  temperature, 
and  flow  steppage  in  the  digester.  All  the  results  are  as 
would  ne  expected.  However  the  working  of  the  model  is  such 
that  to  simulate  a  flow  stoppage  from  the  blow  line  the 
chipmeter  speed  has  to  be  set  to  zero.  The  wash  zone 
results  had  substantial  disagreement  with  those  from 
Johnsson  (49) .  They  state  that  though  the  differences 
between  the  results  are  large,  the  magnitude  cf  error  is  the 
same  for  all  cases.  This  discrepency  is  probably  due  to 
inaccuracies  in  the  mass  transfer  portion  of  the  model. 
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CHAPTER  THREE 


DERIVATION  CF  MODEL 


This  chapter  on  the  Kamyr  digester  model  starts 
with  a  detailed  description  of  the  digester  and  the 
accompanying  equipment  with  typical  operating  conditions. 
This  is  followed  by  the  details  of  a  digester  section  that 
is  representative  of  all  its  zones.  The  assumptions  made 
for  deriving  the  model  and  the  reasons  for  their  choice  are 
discussed  next.  Then  a  review  of  the  kinetic  literature, 
the  derivation  cf  model  kinetics,  and  the  batch  or  cooking 
model  are  presented.  The  chapter  ends  with  reasons  for 
choosing  the  particular  model-  A  comparison  of  Kamyr 
digester  models  and  assumptions  from  Stillman  (93),  Johnsson 
(49)  and  Smith  and  Williams  (88)  and  this  worh  are  given  in 
Table  1- 


3.  1 


Detailed  Digester  Description 


There  are  six  distinct  processing  stages  that  wood 
chips  have  to  pass  through  in  a  continuous  digester  system. 
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namely:  presteaming  and  conditioning,  cooking  liguor 
impregnation,  heating  of  chips  and  liguor,  cooking  of  chips, 
washing  of  chemicals  and  reaction  products,  and  cooling 
before  the  pulp  leaves  the  digester.  A  detailed  piping  and 
instrumentation  diagram  for  the  digester  is  given  in  Figure 
3.  1. 

Chips  are  moved  from  the  chip  bin  to  the  hopper 
that  feeds  the  rotating  chipmeter.  The  chipmeter  rotor 
speed  is  set  to  a  given  digester  wood  feed  rate  and  the  low 
pressure  feeder  then  introduces  the  chips  into  one  end  of 
the  cylindrical  presteaming  vessel.  Here,  they  are  slowly 
conveyed  by  a  rotating  screw  to  the  discharge  end.  The 
vessel  operates  at  15  to  20  psi  and  65  0 C  (150  °F). 

Retention  time  and  temperature  are  limited  to  prevent 
prehydrolysis  of  the  hemicelluloses  for  this  would  limit  tne 
yield  and  the  pulp  beating  and  strength  properties.  The 
presteaming  and  conditioning  stage  thus  accomplishes  removal 
of  non-condensible  gases  and  terpenes,  the  terpenes  are 
separated  later  and  represent  a  valuable  byproduct.  It  also 
raises  the  temperature  and  moisture  content  of  chips  which 
facilitates  liguor  impregnation  and  increases  the  density  so 
that  the  chips  will  sink  more  readily  into  the  digester 


•V 


' 

. 

' 

. 


29 


i"- 

m 


a) 

C 

toO 

> 

cO 

>-1 


a 

<rJ 

M 

Cn 

•H 

Q 

W 

W 

<y 

u 

o 

K 

CU 

U 

Q) 

-»-> 

in 

Q) 

Cr> 

a 

>>1 

s 

rO 

US 


** 


<U 

(-1 

3 

tr> 

•H 

Pm 


30 


cooking  liquor. 

Upon  discharge  from  the  steaming  vessel.  The  chips 
fall  into  the  chip  chute  which  is  partially  filled  with 
liquor  and  impregnation  begins  at  this  time.  The  liquor 
level  of  the  latter  is  held  constant  by  a  controller.  Below 
the  cnip  chute  is  a  rotary  type  high  pressure  feeder  and  the 
chips  are  transfered  from  a  pressure  of  15  psi  to  the 
pressure  at  the  digester  top  of  165  psi.  Chips  and  liquor 
enter  one  of  the  exposed  feeder  pockets  which  are  so  placed 
that  there  is  at  least  one  pocket  in  receiving  and  another 
in  discharging  position.  The  chips  in  the  pocket  are 
retained  while  the  liquor  drains  to  the  intake  of  the  chip 
chute  circulation  pump.  The  chips  are  flushed  out  as  the 
pocket  gets  to  a  horizontal  position  by  liquor  from  the  high 
pressure  top  circulation  pump  and  conveyed  to  the  top 
separator  in  the  digester.  The  top  separator  is  a  helical 
screw  surrounded  by  a  cylindrical  screen  that  separates  the 
liquor,  which  is  returned  to  the  top  separator  circulation 
pump  intake  as  the  chips  pass  into  the  digester  impregnation 
zone.  The  nakeup  cooking  liquor  is  introduced  neiow  the  top 
separator  and  consists  of  white  and  black  liquors  in  the 
correct  alkali  to  wood  (A/W)  and  liquor  to  wood  (1/W) 
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ratios . 

In  the  impregnation  zone  the  chips  rapidly  absorb 
liguor  as  they  move  downwards  by  gravity.  It  is  important 
that  the  chips  be  throughly  impregnated  for  if  impregnation 
is  incomplete,  the  chips  would  be  left  with  uncooked 
centers.  Normally  the  chips  spend  45  minutes  at 
temperatures  ranging  from  100  to  125  °C  (215  to  260  °F) . 

The  cooking  reaction  is  negligible  at  these  temperatures  and 
the  sole  purpose  is  diffusion  of  active  alkali  into  the 
chips. 

Ihe  chips  are  now  ready  for  cooking  and  their 
temperature  is  raised  to  150  °C  (305  °F)  in  the  upper  heater 
zone.  Hot  liguor  enters  the  center  of  the  digester  from  the 
outermost  of  the  concentric  pipes  and  moves  radially 
outwards  to  the  screens  in  the  periphery  of  the  digester- 
Ihere  are  two  screens  but  only  one  of  them  is  extracting  at 
a  time.  This  is  to  allow  the  descending  chip  mass  to  remove 
the  chips  that  are  adhering  to  the  screen  slots  from  the 
preceding  operating  period.  The  switching  of  screens  is 
accomplished  by  a  time  switch  that  opens  control  valves  on 
the  liguor  lines  alternatively-  The  same  operation  is 
repeated  lower  down  in  the  lower  heating  zone  and  the 
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temperature  is  raised  from  150  °C  to  165  °C  (330  °F). 

The  chips  next  move  through  the  cooking  zone  where 
actual  cooking  c£  chips  takes  place.  The  active  chemicals 
in  the  cooxing  liquor,  sodium  hydroxide  (NaCH)  and  sodium 
sulphide  (Na2S)  react  with  the  lignin  and  convert  it  into 
compounds  which  dissolve  in  the  cooking  liguor.  The  liquor 
also  attacks  the  carbohydrates,  which  is  undesirable  for 
this  reduces  the  yield  and  alters  the  pulp  characteristics . 
Cooking  time  is  therefore  selected  to  obtain  the  desired 
yield,  pulp  qualities  and  production  rate.  A  typical  time 
is  about  90  minutes. 

During  passage  through  the  cook  zone  the  pulp 
temperature  rises  by  up  to  ten  degrees  as  the  pulping 
reactions  are  exothermic.  The  pulp  is  then  guenched  by 
injecting  cool  liquor.  The  temperature  of  quench  liquor  has 
a  marked  effect  on  the  reaction  rate  and  if  the  pulp  quality 
is  not  suitable  then  the  cooking  time  is  extended  by 
allowing  cooking  to  continue  into  the  wash  zone.  At  the 
third  set  of  screens,  the  extraction  screens,  the  quench 
liquor  displaces  the  spent  cooking  liquor  through  the  upper 
screen.  The  spent  liquor  flows  to  the  flash  tanks  where 
steam  is  produced  for  the  steaming  vessel.  Liquor  is  then 
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sent  to  the  evaporators  for  recovery-  The  lower  screen 
extracts  dilute  black  liquor  from  the  washing  zone,  in  the 
lower  part  of  the  digester.  The  majority  of  this  liquor  is 
used  for  quenching  and  the  rest  is  mixed  with  the  extraction 
flow  from  the  top  screen  and  passes  to  recovery.  The  quench 
temperature  is  sometimes  independently  adjusted  by  mixing 
with  a  cold  liquor  stream  termed  "cheater”  flow. 

The  chips  continue  downward  and  the  countercurrent 
flow  of  wash  liquor  removes  chemicals  and  dissolved  solids 
from  the  chips  by  displacement  and  diffusion.  This  liquor 
is  introduced  at  the  bottom  of  the  digester  and  is  heated  to 
about  135  °C  (275  °F)  in  the  wash  heater  near  the  bottom. 
This  is  done  to  improve  the  rate  of  diffusion  and  washing 
efficiency  in  the  zone  and  is  termed  "Hi-heat"  washing,  as 
the  temperature  is  higher  than  that  of  the  vacuum  drum 
washers.  Residence  time  in  the  zone  varies  from  60  to  240 
minutes  depending  on  mill  requirements  and  downstream 
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liguor  in  the  preceding  zone.  The  cooled  chips  are  then 
removed  from  the  digester  by  the  outlet  device.  The  chips 
have  to  be  cooled  for  if  removed  at  the  cooking  temperature 
sudden  expansion  and  mechanical  action  results  in  less  of 
pulp  strength.  The  cutlet  device  rotates  at  a  speed  of  3  to 
16  rpm.  The  speed  is  adjusted  to  meet  internal  digester 
conditions  such  as  the  blowing  rate  and  desired  pulp 
consistency.  The  chips  at  a  pressure  of  200  psi  are  removed 
in  three  stages  of  pressure  drop  through  the  blow  unit  to 
the  blow  tank  or  storage. 

3- 2  Cescription  of  a  Digester  Section 

fox  purposes  of  modelling,  the  digester’s  seven 
zones  are  split  up  into  ideal  stirred  tanks  or  lumped 
parameter  sections  of  uniform  temperature  and  composition. 
These  sections  are  not  necessarily  of  egual  volume  as  they 
are  divided  to  match  the  digester  configuration  and  for 
computational  convenience. 

A  digester  section  is  further  subdivided  into 
sub-sections  to  represent  the  three  phases,  oven  dry  wood, 
the  entrapped  or  absorbed  liguor  in  the  wood  and  the  free 
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liquor  flowing  in  the  void  volume.  The  wood  and  the 
entrapped  liquor  are  in  temperature  equilibrium  but  the  free 
liquor  has  a  temperature  of  its  own,  as  external  flews  to 
the  digester  mix  only  with  the  free  liquor,  in  the  wash 
zone  this  liquor  flows  countercurrent  to  the  chip  mass. 

For  material  and  energy  balances  thirteen  state 
variables  are  considered:  Wood  flow  in  mass  of  oven  dry 
wood/minute;  four  wood  components,  lignin,  cellulose,  mannan 
and  xylan  as  mass  fraction  percentages  of  oven  dry  wood;  the 
temperatures  of  the  wood  and  free  liquor;  and  three 
concentrations  each  for  the  entrapped  and  free  liquors, 
sodium  hydroxide,  sodium  sulphide  and  reaction  products. 

3.3  Model  Assumptions 

The  assumptions  made  for  deriving  the  model  are 
presented  and  their  reasons  are  given. 

3-3.1  Steaming  Vessel 

i)  The  temperature  of  chips  and  the  enclosed 
moisture  is  raised  to  that  of  saturated  steam  in  the 
steaming  vessel.  This  is  based  on  condensation  of  steam 
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with  no  steam  losses. 

Ma  thematically , 

Kcond  *Hcond  =  (Ww  *  Cpw  +  Wmo  *  Cpm)  *  (Tsv  -  Tw) 

where 

Wcond  is  the  mass  flow  of  steaming  vessel  condensate, 

Hcond  is  the  latent  heat  of  condensing  steam, 

Kw  is  the  mass  flow  rate  for  wood  chips, 

Wmo  is  the  mass  flow  of  moisture  in  the  wood  chips, 

Cpw  and  Cpm  are  the  specific  heat  of  the  wood  and  moisture, 
Isv  rs  the  steaming  vessel  temperature 
and  Tw  is  the  temperature  of  the  chips. 

ii)  Moisture  in  excess  of  the  fibre  saturation 
point  freezes  when  temperatures  are  below  the  freezing 
point.  The  fibre  saturation  point  is  reached  at  about  30% 
moisture  content. 

The  first  assumption  is  made  so  that  the  steam 
condensate  that  flows  along  with  the  chips  and  dilutes  the 
liquor  can  be  estimated.  This  is  reasonable  as  past 
experience  in  industry  has  shown  (27) .  The  second 
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assumption  has  teen  found  experimentally  by  Stone  and 
Scallan  (95) . 

3-3.2  A  Digester  Section 

i)  Each  digester  sub-section  is  assumed  to  be  an 
ideal  stirred  tank,  i.e.  it  is  perfectly  mixed. 

ii)  The  wood  is  assumed  to  consist  of  lignin, 
cellulose,  mannan  and  xylan-  There  are  also  some  easily 
removable  components,  such  as,  uronic  acid,  acetyl  groups 
and  extractives. 

iii)  The  none  dry  wood  substance  density  is 
constant  at  1-53  (95.5  lb/  cu  ft  )  and  there  is  little 
variation  with  species,  composition  or  type  of  growth. 

iv)  The  chip  movement  is  by  plug  flow  and  the 
velocity  is  a  known  function  of  distance  from  the  digester 
top,  geometry,  pulp  density  and  withdrawal  rate. 

v)  The  chips  undergo  compaction  as  they  move  down 
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the  digester ,  i.e.,  their  bulk  density  increases. 

vi)  Entrapped  liquor  volume  is  based  on  the 
critical  moisture  content  of  the  chips. 

vii)  Digester  liquor  consists  of  sodium  hydroxide, 
sodium  sulphide  and  reaction  products  termed  ’dissolved 
solids  * . 


viii)  The  downward  movement  of  liquor  in  the 
cooking  sections  is  plug  flow  and  the  residence  time  of  both 
the  entrapped  and  free  liquors  is  assumed  to  be  the  same  as 
the  wood. 


ix)  No  back  mixing  of  materials  between  sections  is 

assumed. 


x)  Sodium  sulphide  hydrolysis  is  based  on  an 
empirical  relationship  between  the  ligucr  component 
concentrations  and  the  temperature. 


xi)  All  the  reaction  products  are  converted  to 
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dissolved  solids  and  the  pulp  yield  is  a  sum  of  the  wood 
component  concentrations. 

xii)  Heat  exchange  with  the  digester  walls  is 
divided  proportionately  between  the  wood  and  entrapped 
liquor  and  the  free  liquor. 

xiii)  Entrapped  and  free  liquor  flew  in  their 
respective  sections  with  their  own  temperature  and  heat 
transfer  between  them  depends  on  their  temperature 
difference. 

xiv)  The  liguor  and  wood  specific  heats  and  heat 
transfer  coefficients  are  independent  of  temperature  as 
average  values  in  the  operating  temperature  range  are  used. 

xv)  The  mass  transfer  coefficients  are  equal  for 
the  three  liquor  components. 

xvi)  There  is  no  mass  transfer  of  wood  components 
by  diffusion. 

The  continuous  digester  is  a  tubular  reactor  which 
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can  be  represented  .by  N  stirred  tanks,  especially  when  N  is 
large.  So  assumption  (i)  is  reasonable.  The  basis  for  (ii) 
is  that  the  major  wood  constituents  are  lignin,  mannan, 
xylan  and  cellulose  and  the  easily  removable  components 
dissolve  in  the  early  part  of  the  cook.  The  density  of  wood 
substance,  assumption  (iii) ,  is  from  actual  experimental 
data  by  Stamm  and  Hansen  (91). 

The  residence  time  data  { ( i v) )  was  presented  by 
Hamilton  (34)  and  Johnsson  (49)  from  radioactive  tracer 
studies  on  a  continuous  digester.  The  chip  velocity  is 
observed  to  decreases  down  the  digester  and  this  is  due  to 
chip  compaction  ( (v)  )  . 

The  volume  of  the  entrapped  liquor  ((vi))  in  the 
chips  is  a  function  of  liquor  density,  as  the  chips  swell  up 
and  reach  the  critical  moisture  content  limit  which  is 
known. 

Digester  liquor  composition  {(vii))  is  mainly 
sodium  hydroxide,  sodium  sulphide  and  reaction  products, 
although  chemicals  such  as  carbonates,  sulphates  and  other 
sulphides  are  formed  in  the  recovery  cycle  they  are  either 
negligible  or  inert  to  the  pulping  reaction. 

Plug  flow  for  liquor  (  (viii) )  and  negligible  radial 
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gradients  or  backmixing  ((ix))  are 
mathematical  model.  There  is  also 
chips  (34). 


assumed  to  simplify  the 
no  backmixing  of  the 


(72) 

(93) 


The  sodium 
and  the  actual 


sulphide  hydrolysis  (  ( x ) )  is  from  Martin 
relation  used  is  taken  from  Stillman 


Assumption  ( (xi) )  on  yield  and  reaction  products 
follow  from  a  mass  balance  on  wood.  Since  the  wood 
concentrations  are  mass  fractions  on  a  percent  oven  dry  wood 
basis,  a  sum  of  the  individual  concentra  tions  would  give 
the  percentage  yield. 

Assumptions  (  (xii) -  (xiii) )  are  made  to  simplify  the 
heat  exchange  model  between  the  wood,  free  and  the  entrapped 
liquors.  The  temperature  difference  between  the  wood  plus 
entrapped  liquor  and  the  free  liquor  is  due  to  the  effect  of 
an  external  heater  flow,  which  mixes  with  only  the  free 
liquor.  It  takes  the  wood  and  entrapped  liguor  some  time  to 
reach  that  temperature. 

The  liguor  and  wood  heat  transfer  coefficients  are 
taken  to  be  constant  (  (xiv) ) ,  as  data  was  not  available  for 
their  variation  with  temperature.  The  wood  specific  heat 
relationship  with  temperature  has  been  studied  by  Dunlap 
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(22)  and  an  average  value  in  the  digester  operating  range  is 
used  in  the  model-  The  liquor  specific  heat  value  was 
similarly  determined.  The  assumption  for  diffusion  of  wood 
component  (  (xvi) )  being  neglected  is  because  the  transport 
of  pure  substances  by  the  flow  of  wood  and  liquor  is  much 
larger  than  that  by  diffusion. 


3. 3.  3 

sect  ion 
leaving 


Impregnation  Zone 

i)  The  flow  of  chips  entering  the  first  digester 
is  based  on  the  chip  meter  flow  and  the  dIow  flow 
the  digester,  and 


ii)  The  easily  removable  wood  components  are 
assumed  to  completely  react  in  the  impregnation  zone. 

The  first  assumption  {  (i) )  for  chip  flow  to  the 
first  section  is  apparent  when  a  stoppage  in  chipmeter  or 
blow  flow  is  considered. 

The  easily  removable  wood  components  have  a  high 
reaction  rate  and  it  is  observed  that  they  are  not  present 
at  the  end  of  the  cook.  Simulation  studies  by  Stillman 
(93),  who  considered  no  such  assumption  (  (ii)  )  showed  this 


to  be  true 
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3.3.4  Heating  Zone 

i)  Heater  flows  involve  no  mass  transfer  and  only 
help  to  raise  the  free  liquor  temperature. 

This  assumption  is  to  simplify  calculations,  since 
the  external  flow  is  a  pump  around  stream  and  would  have  the 
component  concentration  equal  to  that  in  the  digester 
section  from  which  it  is  removed,  at  steady  state. 

3.3.5  Cooking  Zone 

i)  Lignin,  mannan  and  xylan  have  varying  reaction 
rates,  depending  on  their  concentration  and  process 
conditions.  These  rates  are  termed  extractive,  bulk  and 
residual  and  relate  to  the  ease  of  removal  of  the  wood 
components. 

ii)  Tne  reaction  rate  for  a  wood  component  Z  is  of  the  first 

order  type  and  can  be  written  as: 

dZn  =  -kij*Zn 
dt 

where  Zn  is  the  concentration  of  component  Z  in  section  n 
and  kij  is  the  reaction  velocity  constant  that  is  dependent 
on  the  liquor  concentrations  and  temperature. 
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iii)  The  reaction  rate  is  dependent  on  the  entrapped  liquor 
concentration. 

Assumptions  (i)  to  (iii)  and  comments  on  them  are 
presented  in  section  3, 4  , 

3-3-6  Extraction  /  Quench  Zones 

The  extraction/quench  screens  with  their 
complicated  flew  patterns  are  treated  as  a  perfect  mixing 
zone  for  free  liquor,  from  the  cooking  and  wash  zones-  The 
flow  patterns  assumed  are  shown  below. 
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This  assumption  for  perfectly  mixed  free  liquor 
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streams  fellows  from  the  earlier  assumption  for  heater  flows 
(3.3.2  (i) ) . 

3.3.7  Washing  Zone 

Free  liquor  in  the  wash  zone  is  countercurrent  to 
the  wood  flow  but  there  is  considered  to  be  no  backmixing 
between  sections.  Plug  flow  is  assumed. 

This  assumption  is  just  an  extension  of  3.3.2 
(viii)  r  3.3.2  (ix)  and  3.3.2  (xvi)  . 


Table  1  Comparison  of  Kamyr  Digester  Models 
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3-4  Model  Kinetics 

A  literature  survey  and  review  of  past  kinetic 
models  is  presented,  followed  by  the  kinetic  model  used  in 
this  work  and  the  reasons  for  its  choice. 

3.4.1  Literature  Survey 

The  pulping  reaction  and  its  kinetics  are  very 
complex  and  have  been  studied  by  a  number  of  researchers  but 
there  is  still  no  general  consensus  on  a  single  uniform 
model  (59,  67,  101).  Pulping  is  governed  mainly  by  the  rate 
of  delignif ication  and  this  reaction  has  been  extensively 
examined.  The  undesirable  side  reactions,  such  as 
carbohydrate  degradation  cannot  be  avoided  and  have  to  be 
minimized  sc  their  kinetics  have  been  of  interest  also. 

Lignin  occurs  in  amounts  varying  from  20  to  35%  of 
the  wood  depending  on  the  species.  It  is  a  large  three 
dimensional  polymer  made  up  of  phenyl-propane  monomers  that 
are  linked  together  in  various  forms.  A  suggested  reaction 
mechanism  is  at  the  solid-liquid  interface  by  the  hydroxide 
or  sulphide  (or  hydrosulphide)  ions  present  in  the  liquor 
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(59) *  These  fragments  dissolve  as  phenolate  or  carboxylate 
anions  or  they  can  undergo  reactions  with  themselves  ,  the 
undissolved  lignin  ,  or  possibly  also  with  the 
carbohydrates.  The  lignin  that  undergoes  alkaline 
condensation  reactions  is  more  difficult  to  remove  and  so 
this  reaction  has  to  be  prevented.  There  are  indications 
that  sulphide  or  hydrosulphide  ions  help  in  reducing  these 
condensation  reactions  (59) . 

In  most  alkaline  pulping  the  desired  speci f ications 
are  met  by  50%  lignin  removal,  but  along  with  this  10  -  15% 
of  the  cellulose  and  at  least  50%  he nicellulose  are  removed. 
As  mentioned  earlier,  the  cellulose  removal  is  undesirable 
for  it  results  in  approximately  a  55%  pulp  yield  based  on  an 
oven  dry  wood  basis  and  is  a  very  important  economic  factor 
as  the  cost  of  wood  is  a  major  part  of  the  total  production 
cost  of  pulp.  For  this  reason  efforts  have  been  made  to 
minimize  the  carbohydrate  reactir  and  enhance  the  yield. 

In  the  early  stages  of  the  cook,  ^emice llulcses  are  reduced 
by  40%  while  the  lignin  dissolved  is  just  20%  (59) .  The 
introduction  of  certain  additives  such  as  anthraguinone  to 
the  pulping  liguor  look  promising  in  reducing  the 
carbohydrate  reaction.  Recently  patents  have  been  issued 
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for  methods  of  reprecipitating  the  hemicelluloses  that  is 
dissolved  in  the  earlier  stages  of  the  cock  (46)  . 

The  carbohydrates  consist  of  basic  linear  polymers 
made  up  of  glucose  (cellulose) ,  xylose  (xylans)  and  a 
copolymer  of  glucose  and  mannose  (glucomannans  and 
galactoglucomannans) .  The  carbohydrate  polymer  degradation 
process  is  dependent  on  polymer  size  or  degree  of 
polymerization  and  the  accessibility  to  pulping  liquor. 

There  are  predominiantly  three  basic  types  of  reactions. 

The  first  is  the  “peeling"  reaction,  which  causes  the  bulk 
of  the  yield  loss,  and  is  the  removal  of  the  reducing  end 
groups  from  the  polymer.  After  an  end  group  has  been 
removed,  the  remaining  end  group  can  again  be  removed  by 
"peeling"  or  can  be  stabilized  by  conversion  to  an  acid 
group.  This  is  the  second  type  and  is  referred  to  as  the 
"stopping"  reaction.  Chemically  the  "peeling"  reaction 
relates  to  the  free  carbonyl  groups  which  are  not  stable  in 
alkali  and  split  up  with  about  sixty  monomer  end  groups  as 
iso-sacchar inic  acids.  In  the  "stopping"  reaction  the 
molecule  is  stabilized  by  formation  of  meta-saccnar inic  acid 
end  groups.  The  third  reaction  class  is  "random  cleavage" 
of  the  molecule.  This  occurs  at  temperatures  of  above  150  C 


- 


54 


with  random  alkaline  hydrolysis  of  accessible  gluco 
linkages.  This  exerts  a  dual  effect  in  the  loss  of 
for  small  alkali  soluble  polymers  are  created  and  a 
carbonyl  end  groups  are  formed  where  further  "peeli 
reactions  can  occur. 

Overall,  glucomannans  and  hardwood  xylans 
especially  susceptible  to  the  ’'peeling"  reaction, 
former  is  the  main  hemicellulosic  component  in  soft 
the  latter  in  hardwoods.  Softwood  xylans  are  relat 
stable  and  only  25  -  30%  of  them  are  removed  in  nor 
pulping  (59).  Several  theories  exist  but  it  is  not 
why  this  is  so  (70  p431#  81  p604) . 

The  next  important  pulping  variable  is  tne 
composition.  Normal  kraft  liquor  is  made  up  of  an 
solution  of  sodium  hydroxide  and  sodium  sulphide  wh 
sodium  sulphide  is  introduced  as  sodium  sulphate  in 
chemical  recovery  cycle.  This  is  due  to  the  faster 
reaction  rate  observed  with  the  addition  cf  sodium 
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Ihe  kraft  liquor  components  have  the  following  equilibria 
(101)  : 


NaOH 

(Na+)  + 

(OH") 

(i) 

Na  S 

2 

2  (Na  "*)  + 

(S2-  ) 

(ii) 

(S2-) 

+  (OH")  + 

(HS“) 

(iii) 

NaSH 

r (Na+ )  + 

(HS~) 

(iv) 

(HS“) 

+  H  0  k - (OH") 

♦  H  S 

2 

(v) 

Liquor  concentration  definitions  involve  these 
equilibrium  concepts.  The  active  alkali  definition 
assumes  reactions  (iii)  and  (v)  are  complete  which  means 
that  the  hydroxide  ion  concentration  is  a  sum  of  the  sodium 
hydroxide  and  sodium  sulphide  molar  concentrations.  While 
the  effective  alkali  definition  assumes  that  reaction  (iii) 
is  complete  and  that  reaction  (v)  is  insignificant.  The 
alkali  concentration  is  then  a  sum  of  the  sodium  hydroxide 
and  half  that  of  the  sodium  sulphide. 

It  is  difficult  to  determine  the  exact  hydroxide, 
sulphide  or  bisulphide  equilibria  involved  especially  the 
effect  of  the  individual  contribution  of  sulphur  ions 
towards  pulping  (101) ,  but  consistent  pulping  results  are 


- 

. 


56 


obtained  if  the  effective  alkali  concept  is  used  as  opposed 
to  that  of  the  active  alkali. 

3.4.2  Review  of  Kinetic  Models 


The  overall  pulping  reaction  involves  five  steps 
which  proceed  in  sequence-  If  one  step  is  slew  compared  to 
the  others  then  it  is  said  to  control  the  overall  rate. 

These  steps,  with  respect  to  pulping  are  movement  of 
reactants  or  chemicals  into  the  chips,  adsorption  of 
chemicals  onto  active  sites,  the  actual  chemical  reaction, 
desorption  of  reaction  products  and  transport  of  desorbed 
reaction  products  out  of  the  chip.  The  first  and  the  fifth 
transport  steps  involving  the  transport  of  chemicals  and 
reaction  products  become  more  prominent  as  the  distance  from 
the  chip  surface  increases.  This  was  observed  from  studies 
on  chip  thickness  and  shape  on  pulping  rate.  Backman  (3) , 
Hartler  and  Ostberg  (35) ,  Larocgue  and  Maass  (63)  and 
Kulkarni  and  Nolan  (61)  conducted  such  experiments  and  found 
that  the  thickness  was  important  above  a  certain  "critical 
value"  which  varied  with  the  cooking  temperature-  This 


critical  thickness  has  been  assigned  varying  values  of  1  mm 
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(3),  3  mm  (35)  and  up  to  10  mm  (63)  but  the  effects  of 
diffusion  can  be  minimised  if  a  chip  with  a  thickness  of  2 
to  3  mm  oi  less  is  used  with  current  digester  pulping 
temperatures  (101).  The  same  observations  are  obtained  when 
the  reaction  rate  coefficient  of  30  Kcals/mcle  (21)  is 
compared  with  that  of  diffusion  limited  processes  at  5  -  6 
Kcal/mole  (81  p3G9)  .  For  soda  pulping,  Daleski  (21) 
suggests  195  °C  as  the  temperature  above  which  diffusion 
becomes  dominant,  liquor  diffusion  is  important  in 
impregnation  (81  p265)  and  additional  methods  of  forced 
penetration  by  liquor  and  gas  flow  or  high  pressure 
impregnation  have  been  developed  to  achieve  uniform 
impregnation  and  thereby  uniform  pulping  with  minimum  amount 
pulp  rejects.  The  effects  of  diffusion  on  pulping  were  also 
studied  by  Jchnsson  (49)  using  computer  simulation  and 
varying  chip  thickness  distribution,  penetration  depth  and 
diffusion  coefficients,  and  similar  results  on  pulping  and 
rejects  were  obtained- 

The  chemical  reaction  step  involves  two  dissimilar 
phases,  wood  and  liquor  which  gives  a  heterogeneous  reaction 
system.  This  is  extremely  difficult  to  evaluate  and  so 
homogeneous  models  have  been  generally  derived 
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mathematically  in  terms  of  a  rate  expression  and  from 
analysis  of  actual  cooking  data.  Several  models  have  been 
developed  on  this  basis  and  the  effects  of  liquor 
composition  and  temperature  included.  Most  of  the  earlier 
work  was  done  for  soda  liquor  which  is  a  limiting  case  of 
kraft  liquor,  since  only  sodium  hydroxide  is  present  as  a 
pulping  agent.  The  soda  process  was  used  extensively  to 
pulp  hardwoods  which  are  easier  to  pulp  and  hence  require 
less  severe  pulping  conditions.  This  is  due  to  a  lower 
lignin  content  in  hardwoods.  The  introduction  of  the  kraft 
process  gave  a  higher  pulping  rate  and  softwoods  too  could 
then  be  pulped  economically.  Though  woods  have  been  broadly 
isolated  into  hardwood  and  softwood  classes,  each  individual 
species  has  its  own  reaction  rate  (81  p614). 

A  first  order  kinetic  representation  is  dominant 
and  Mitchell  and  Yorston  (76),  Schwartz  and  Bray  (84), 
Larocque  and  Maass(63)#  Kulkarni  and  Nolan  (61)  and  Carroll 
(14)  have  derived  delignif ication  equations  of  the  form: 

di  =  ki*  (L  -  Ld)  3.4.  (i) 

dt 

where  ki  is  the  kinetic  constant,  L  the  lignin  present 
at  the  time  t  and  Ld  is  lignin  that  is  difficult  to  remove 
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in  the  last  stages  of  the  cook. 

The  kinetic  constant,  k  is  related  to  temperature 
ty  an  Arrhenius  type  relationship  and  activation  energy 
values  vary  from  30,000  to  35,000  cal/mole  for  kraft  cooks- 
The  effect  of  hydroxide  and  sulphide  concentration  in  the 
liquor  is  included  in  the  effective  rate  constant 
calculation  as  a  linear  relationship. 

These  earlier  kinetic  models  related  to  only  the 
"bulk”  reaction.  Kleppe  (59)  and  Stillman  (93)  state  that 
the  pulping  reaction  occurs  with  three  different  rates, 
"extractive",  "bulk"  and  the  "residual",  while  Kieinert  (55, 
58)  mentions  only  the  "bulk"  and  "residual"  rates.  The 
"extractive"  rate  corresponds  to  the  easily  removable  and 
readily  accessible  wood  constituents,  the  "bulk"  to  describe 
the  majority  of  the  reaction  and  the  "residual"  for  the  last 
few  percent  of  least  accessible  and  most  difficult  amount  of 
wood  constituent  to  remove.  These  differing  rates  are 
assumed  to  be  first  order.  Stillman  derived  a  relation  to 
identify  the  breakpoint  between  bulk  and  residual  lignin  as 
a  function  of  lignin  and  liquor  concentrations  and  the 
cooking  temperature.  The  extractive  rate  is  based  on 
extrapolating  the  bulk  reaction  curve  to  zero  cooking  time 


* 


. 


. 


60 


and  noting  the  difference  between  the  actual  and  bulk  lignin 
content  obtained  by  extrapolating  the  curve.  Stillman 
plotted  data  from  Hobden  (41),  Kleppe  (59),  and  Aurell  and 
Hartler  (1,  2)  to  show  the  three  rates. 

Mitchell  and  Yorston  (76)  ,  Stillman  (93)  and  Wilder 
and  Daleski  (102)  too  used  a  first  order  equation  of  the 
form  : 

dl  -  -  k*L  3.4. (ii) 

dt 

to  describe  the  reaction  kinetics.  The  rate  constant  k  is  a 

function  of  both  the  kraft  liquor  component  concentrations 

k  =  k  *  (OH  ")  +  k  *(S2")**a  3.4.  (iii) 

1  2 

where  (OH")  is  the  effective  sodium  hydroxide  concentration, 

(S2")  is  the  sulphide  concentration  and  ’a*  is  the 

exponential  power  for  sulphide  which  varies  with  wood 

species.  Daleski  (21)  obtained  38  Kcal/mole  for  and  26.9 

Kcal/mole  for  k  . 

2 

Stillman  uses  the  same  form  of  equation  for  pulping 
of  seven  hardwood  and  softwood  components.  The  effective 
alkali  relationship  is  based  on  data  from  Martin  (72). 

Since  sodium  sulphide  was  found  to  have  a  large  effect  on 
carbohydrate  removal  in  the  initial  portions  of  the  cook,  a 
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effective  concentration  for  sodium  sulphide  is  evaluated  to 
obtain  a  good  correlation.  To  calculate  the  three  different 
reaction  rates  for  bulk,  extractive  and  residual  phases, 
factors  which  multiply  the  Arrhenius  constants  are 
calculated. 

.Findley  and  Nolan  (25)  derived  a  model  for  two 

simultaneous  reactions  for  wood  components  in  neutral 

sulphite  pulping.  They  are: 

dYa  =  -Ka*Ya 
d  t 

dYb  =  -Kb* Yb  3-4.  (iv) 

dt 

and  Y  =  Ya  *  Yb  +  43 

where  Ya  and  Yo  are  wood  fractions  with  velocity  constants 
Ka  and  Kb  respectively  and  Y  is  the  wood  yield  with  43% 
unreacted  wood.  Johnsson  (49)  uses  a  somewhat  similar  model 
for  his  kinetics. 

The  latest  models  are  from  LeMcn  and  Teder  (67)  and 
Smith  and  Williams  (88)  and  they  use  the  following  equation 


-  dL  =  K  *L*  (0 H")  +  K  *1*(0H")*(HS  )  **b  3.  4.  (v) 

dt  1  2 

where  (OH-)  is  the  effective  alkali  and  (HS  )  the  bisulphide 
concentration.  The  exponential  constants  used  by  Smith  and 
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Williams  are  unity  and  they  have  used  the  same  form  of 
equation  for  five  wood  components:  fast  and  slow  reacting 
lignin,  cellulose,  xylan  and  mannan.  Kinetic  data  is 
derived  from  actual  laboratory  cooks  of  hardwood  and 
softwood  species.  LeMon  and  Teder  have  used  the  equation 
only  for  lignin  and  have  fractional  values  for  the 
exponential  constants  rather  than  unity. 

3.4.3  Kinetic  Model 

The  kinetic  model  in  this  work  is  based  on 
Stillman's  (93)  model.  Therefore  the  kraft  liquor 
components,  sodium  hydroxide  and  sodium  sulphide  are 
considered  as  actual  hydroxide  and  sulphide  in  the  kinetic 
equation  instead  of  effective  alkali  and  bisulphide.  The 
total  hydroxide  or  alkali  concentration  is  taken  as  the 
sodium  hydroxide  in  the  liquor  plus  the  amount  generated 
from  the  sodium  sulphide  hydrolysed.  For  the  sodium 
sulphide,  instead  of  assuming  fifty  percent  hydrolysis  as  is 
defined  for  the  effective  alkali,  Martin's  (72)  data  for 
hydrolysis  is  used  and  the  actual  hydrolysis  calculated. 
Sodium  sulphide  was  considered  as  the  sulphide  although 
hydrolysis  generates  bisulphide  or  hydrogen  sulphide.  An 
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effective  sodium  sulphide  concentration  is  evaluated  fox 
degradation  of  carbohydrates,  as  the  sulphide  is  responsiole 
for  a  higher  pulping  reaction  in  the  earlier  parts  of  the 
cook  (59,  101). 

The  reaction  kinetics  are  considered  first  order 
with  respect  to  the  wood  components  and  the  contribution  of 
sodium  hydroxide  and  sodium  sulphide  to  the  kinetic  constant 
additive.  The  relevant  equations  follow: 

dXn  =  -Ko*Xn  3-4.  (vi) 

dt 

where  Ko  is  the  reaction  constant  and  Xn  the  wood  component 
concentration  in  the  nth  digester  section.  The  constant  Ko 
is  equal  to: 

Ko  =  K  *  (OH  “)  +  K  *(s2-)**a  3-4.  (vii) 

1  2 

The  exponent  *a'  is  0.631  for  hardwoods  and  0.681  for 
softwoods  (93).  The  constants  K ^  and  K2  are  for  the 
hydroxide  and  sulphide  respectively  and  related  to  the 
cooking  temperature,  T,  by  the  Arrhenius  equation  and  a 

reaction  phase  factor,  B. 

eg.  K  =  B  *exp  (A  -  E  /T)  3.4  (viii) 

11  11 

The  constants  A,  B  and  E  for  various  wood  components  are 


shown  in  Table  2. 
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Lignin  and  xylan  are  assumed  to  react  with  three 
different  rates,  extractive,  bulk  and  residual  while  mannan 
is  considered  to  have  only  extractive  and  bulk  rates-  The 
extractive  amount  is  calculated  by  extrapolating  the  bulk 
reaction  curve  to  zero  cook  time  and  the  difference  between 
the  component  concentration  thus  obtained  is  the  actual 
value  assumed  to  re  in  the  fast  extractive  phase.  The 
breakpoint  between  the  bulk  and  the  residual  phases  is 
calculated  from  the  wood  compositions  and  the  cooking 
conditions,  liquor  concentrations  and  the  temperature  (93) . 
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3.  4.4 


Reasons  for  Kinetic  Model 


The  reaction  model  described  above  was  chosen 
because  it  seems  to  match  the  physical  realities  of  the 
cooking  process.  The  complex  sulphur  ion  hydrolysis  was  not 
used  as  it  is  difficult  to  determine  the  exact  guantities 
involved.  A  fifty  percent  sulphide  hydrolysis  concept  of 
effective  alkali  was  not  used  for  similar  reasons,  instead 
sulphide  hydrolysis  is  calculated  and  the  alkali 
concentration  evaluated  in  terms  of  measurable  liquor 
components.  The  differing  rates  of  reaction  observed  for 
some  wood  components  are  included  as  cooking  does  sometimes 
go  into  the  residual  reaction  phase,  although  this  is  not 
desirable  as  the  yield  loss  at  this  stage  is  greater  than 
the  nenefit  gained  by  the  extra  lignin  removal.  This  has 
been  observed  and  reported  by  several  authors,  as  mentioned 
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Table  2  Kinetic  Rate  Constants  for  Softwood  and  Hardwood  Species 
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CHAPTER  FOUR 


DESCRIPTION  CF  MODEL 

The  Kamyr  continuous  digester  model  derived  in  tnis 
work  could,  based  on  mathematical  modelling  principles,  have 
been  anywhere  between  two  extremes.  A  set  of  trivial  linear 
algebraic  equations  could  have  been  derived  for  the  digester 
variaoles  using  input-output  plant  data.  This  may  have 
accomplished  the  gob  but  only  within  a  certain  limited 
operating  range  and  thus  have  severe  restrictions.  A 
rigorous  model  expressed  in  terms  of  several  nonlinear 
partial  differential  equations,  would  have  been  mainly  a 
mathematical  and  computing  exercise  and  would  have  served 
little  purpose  in  what  this  project  was  eventually  aimed  at, 
i.e.  control  of  the  cooking  process.  An  intermediate 
approach  was  taken  instead  and  this  helped  in  making  the 
mathematical  and  computational  requirements  of  the  model 
reasonable  but  without  compromising  the  physical  realities 
of  the  system. 

In  this  chapter  the  derivation,  solution  technique, 
computer  programming  and  testing  of  the  digester  model  are 
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presented . 

4.1  Derivation  of  Model  Equations 

Since  the  continuous  Kamyr  digester  is  a 
distributed  parameter  system,  a  partial  differential 
equation  was  derived  for  a  general  wood  component 
concentration.  The  equation  is  presented  in  Appendix  A-1- 
This  equation  was  reduced  to  an  ordinary  differential 
equation  which  applies  to  each  of  the  N  sections  that  the 
digester  is  divided  up  into.  This  makes  each  section 
a  ‘lumped  parameter*  or  'stirred  tank*  system.  The 
remaining  equations  are  derived  directly  fcr  such  a  system. 
The  alternate  approach  would  have  been  to  directly  start 
from  the  concept  of  equivalence  of  a  tubular  reactor  to  a 
series  of  stirred  tanks  and  to  develop  the  model  equations 
from  there  onwards.  This  was  not  done  as  the  existing 
mathematical  derivation  makes  the  set  of  boundary  conditions 
for  the  discretized  case  more  obvious. 

Material  and  energy  balances  for  thirteen  state 
variables  were  derived  for  an  ideal  stirred  tank,  giving  a 
set  of  thirteen  coupled,  first  order,  nonlinear  ordinary 
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differential  equations.  This  set  of  equations  is  shown  in 
Table  3  and  again  in  Appendix  A. 2. 

The  partial  differential  equation  (Eqn  A.  1-1)  was 
simplified  ty  neglecting  the  radial  and  angular  components 
and  discretizing  along  the  vertical  axis.  Eiffusion  of  the 
solid  phase  wood  components  is  assumed  negligible,  which 
eliminates  the  second  order  terms,  leaving  convection,  toxal 
reaction  and  unsteady  state  terms.  Discretization  of  the 
initial  conditions  led  to  an  initial  entering  condition  at 
tne  top  of  the  digester  and  an  initial  profile  of  starting 
values  for  each  of  the  N  -digester  sections. 

The  first  eguation  in  the  set  of  equations  is  for 
the  wood  density  in  a  digester  section  and  this  is 
essentially  to  calculate  the  wood  flow  rate  leaving  the 
section.  The  flow  rate  into  the  first  section  depends  on 
the  chipmeter  speed  and  the  blow  flow  out  of  the  digester. 
The  flow  out  of  a  digester  section  depends  on  the  entering 
flow  and  the  wood  density  in  that  section.  The  wood  density 
is  related  to  the  degree  of  compaction  which  is  a  measure  of 
the  reduction  in  wood  volume  per  unit  mass.  The  compaction 
relationship  is  a  dynamic  one  and  is  based  on  the  amount  of 
cooking  that  the  chips  have  been  subjected  to  in  the 
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digester,  i.e.  their  H-factor.  The  next  four 
the  wood  component  balances  for  lignin,  cellulo 
and  xylan.  These  equations  include  flow  by  con 
their  depletion  by  first  order  kinetics  that  we 
the  last  chapter.  Liquor  component  balances  in 
entrapped  liquors  are  for  the  hydroxide,  sulphi 
reaction  products  and  include  convection  flew  a 
transfer  of  the  components  between  the  entrappe 
liquor  volumes  in  the  section.  The  depletion  o 
chemicals  is  calculated  from  the  rate  of  reacti 
constituents  and  their  corresponding  stoichiome 
coefficients  and  is  deducted  from  the  entrapped 
mass  balance  for  the  reaction  products  or  disso 
in  the  entrapped  and  free  liquors  contain  terms 
the  chemical  balances.  The  rate  of  generation 
dissolved  solids  is  equal  to  the  rate  of  deplet 
and  cooking  chemicals. 

The  energy  balances  are  evaluated  for 
entrapped  liquor,  and  the  free  liquor.  Here  mo 
are  made  in  the  free  liquor  calculations  to  acc 
external  heater  and  extraction  flows.  Perfect 
external  stream  and  the  free  liquor  is  assumed. 
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Table  3  Nonlinear  Digester  Model  Eyuationsi 
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convective  transfer  between  the  wood  and  free  liquor  phases 
are  included.  The  heat  of  reaction  generated  and  the  heat 
transfer  from  the  walls  are  divided  proportionately  among 
the  two  equations. 

4.2  Equation  Solution  Technique 

The  digester  model  was  programmed  in  FORTRAN  and 
run  on  an  Amdahl  470  V/6  computer.  Simplified  flow  diagrams 
of  the  program  are  shown  in  Figure  4.1,  and  in  Appendix  C.l. 
Detailed  flow  diagrams  and  program  listings  are  included  in 
a  separate  programmers  manual.  Starting  values  for  an 
initial  profile  are  read  in,  estimated  based  on  the  entering 
and  operating  conditions  or  calculated  as  a  steady  state 
profile.  The  first  half  of  the  program  for  cocurrent  chip 
and  liquor  flow  in  the  impregnation,  heating  and  cooking 
zones  of  the  digester  is  an  initial  value  problem.  The 
second,  or  lower  part  of  the  digester  which  includes  the 
quench,  extraction,  wash  and  cold  blow  zones  have  a  wash 
flow  countercurrent  to  the  chip  flow.  This  problem  can  be 
treated  as  an  initial  value  problem  also,  as  long  as  the 
starting  values  and  the  profile  are  known.  The  steady  state 
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profile  for  the  countercurrent  zones  has  to  be  solved 
iteratively  as  a  boundary  value  problem. 

The  initial  value  problem  was  integrated  using  an 
explicit  fourth  order  Runge-Kutta  technique  (  Runge-Kutta 
Merson  )  with  truncation  error  estimation  and  integration 
step  size  control.  The  subroutine  used  is  DASCRU,  from  tne 
IMSL  (International  Mathematical  and  Statistical  Library, 
(45))  subroutine  package  and  was  found  to  be  the  best 
available.  A  comparison  of  computing  times  and  accuracy  for 
other  solution  methods  is  given  in  Appendix  B.l. 

The  steady  state  digester  profile  which  is  the 
solution  of  a  set  of  nonlinear  simultaneous  algebraic 
equations  is  calculated  using  a  modified  form  of  the 
Newton-Eaphson  method.  This  method  is  called  Brown's  method 
and  is  available  in  the  IMSL  package  as  the  subroutine, 
ZSYSTM  (45).  For  the  cocurrent  problem  a  set  of  thirteen 
equations  are  solved  simultaneously  to  obtain  the  steady 
state  values  of  the  digester  variables  in  ore  section,  but 
in  the  washing  zone,  equations  for  all  the  sections  have  to 
be  setup  and  solved  simultaneously. 
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4.3  Description  of  Computer  Propram 

The  computer  program  for  the  Kamyr  digester  model 
solves  the  set  cf  nonlinear  ordinary  differential  equations 
shown  in  Table  3  for  each  section  of  the  digester.  The 
overall  program  is  a  dynamic  one  and  allows  for  cooking  in 
the  wash  zone. 

The  program  starts  with  real  time  being  initialized 
and  all  the  relevant  flags  set.  The  user  defined  flags 
direct  the  logical  flow  of  the  program  according  to  the 
options  specified.  For  example,  flags  can  be  set  to  specify 
omission  of  heat  balances  {  which  reduces  the  number  of 
equations  tc  be  solved  ),  plotting  of  profiles,  which 
variables  are  to  be  plotted  and  chip  history.  Input  data 
that  is  read  in  the  mainline  program  includes  physical 
constants,  initial  digester  entering  conditions,  digester 
configuration  data  and  an  initial  digester  profile,  if  it  is 
known.  The  zone  and  section  volumes,  internal  piping  and 
digester  wall  areas  are  all  calculated  for  the  specified 
digester  configuration.  These  areas  and  section  volumes  are 
printed  out  along  with  the  section  numoers  which  have 


. 


79 


screens  for  liquor  extraction  in  them.  A  sample  printout  is 
shown  in  Appendix  C.2.  The  initial  digester  profile  is 
calculated  next,  if  it  has  not  been  specified  as  a  part  of 
the  input  data. 

The  main  processing  loop  starts  next  with  two 
objectives:  to  calculate  a  digester  profile  at  a  certain 
time  and  to  keep  track  of  each  entering  wood  mass  so  as  to 
trace  its  cooking  history  through  the  digester. 

The  integration  step  size  is  specified  and  the 
input  conditions  to  a  section,  at  time  t,  are  set  to  the 
digester  entering  conditions  in  the  first  section  or  the 
initial  condition  of  the  previous  section  at  the  same  time. 
The  outlet  conditions  from  the  section  are  set  to  the 
initial  conditions  of  the  section  at  time  t.  The  set  of 
equations  are  then  solved  as  an  initial  value  problem  and 
conditions  at  time  t  +  At  calculated.  The  section 
temperatures  are  estimated  by  assuming  them  to  be  equal  to 
the  ones  in  the  last  time  interval.  Next,  a  check  is  made 
whether  the  section  has  a  screen  and/or  an  external  entering 
flow  such  as  one  of  the  heater  streams,  and  if  so,  a  perfect 
mixing  calculation  is  made  for  the  chemical  concentrations 
and  temperature  of  the  entering  free  liquor  stream  to  the 
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section.  The  rate  constants,  mass  transfer  coefficients, 
compaction  in  the  section,  wood  mass  and  free  and  entrapped 
liquor  volumes  are  calculated  and  the  thirteen  first  order 
differential  equations  set  up.  The  equations  are  solved  b y 
a  Bunge  Kutta  type  routine  from  the  IMS!  subroutine  package. 
If  the  section  temperature  calculated  from  the  energy 
balance  is  within  a  specified  limit  of  the  assumed  value, 
then  the  wood  conditions  are  updated.  Otherwise  the  section 
temperatures  are  iterated  upon  after  recalculating  the 
temperature  dependent  variables.  For  the  countercurrent 
flow  sections,  the  starting  values  are  estimated  first  and 
only  the  steady  state  problem  is  solved  differently.  The 
dynamic  boundary  value  problem  is  converted  to  an  initial 
value  problem  using  the  known  countercurrent  stream  entering 
conditions  or  the  calculated  'steady  state1  profile  whicn 
has  been  solved  as  a  boundary  value  problem. 

After  calculations  are  made  for  all  digester 
sections,  the  profile  is  printed  out  and  before  calculations 
are  repeated  from  the  top,  the  current  entering  or  process 
operating  conditions  are  checked.  If  they  are  different,  a 
new  set  of  conditions  are  calculated. 

At  this  stage  further  calculations  and  execution 
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are  stopped  if  one  of  the  stopping  criteria  are  met  which 
are  time,  number  cf  profiles  or  when  the  character  string 
*  STOP  *  is  read  in. 

4.4  Testing  of  Model 

The  model  was  tested  in  three  parts.  First,  the 
working  of  individual  subroutines,  next  numerical  comparison 
of  integration  methods  and  last  the  testing  of  the  overall 
program. 

The  program  subroutines  were  individually  tested 
against  published  data  and  found  accurate  for  modelling 
purposes  within  the  digester  operating  range.  Routines 
tested  first  were  those  that  could  be  checked  with  a  hand 
calculator.  Such  routines  were  the  ones  for  the  digester 
configuration,  residence  time,  digester  entering  conditions 
and  external  flow  calculations.  The  water  density 
calculation  routine  which  is  used  to  evaluate  liguor  density 
change  versus  temperature  was  checked  in  a  temperature  range 
of  25  to  200  °C .  The  results  were  as  expected  since  a 
second  order  polynomial  was  fitted  using  data  values  in  the 
same  range.  The  cooking  liquor  contains  sodium  sulphide 
that  hydrolyses  to  sodium  hydroxide.  The  results  from  the 
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hydrolysis  routine  were  compared  with  data  from  Martin  (72)  - 
These  results  were  reasonable  at  most  points  but  had  a 
maximum  error  of  12%.  The  largest  errors  were  at 
temperatures  around  25  °C  and  at  high  temperatures  and 
sulphidities  which  are  out  of  the  usual  operating  range. 

The  integration  problem  for  the  differential 
equations  was  such  that  a  very  efficient  routine  was 
required.  Fourth  order  Runge  Kutta  and  Predictor  Corrector 
methods  were  examined.  These  routines  were  taken  from  the 
IMSL  (45) ,  IEM  Scientific  Subroutine  Package  (83)  and  some 
were  written  as  a  part  of  the  whole  digester  modelling 
package.  Computer  run  times  and  accuracy  using  the  IMSL 
routine  were  not  as  good  as  with  the  fourth  order  Runge 
Kutta  routine,  when  tested  against  simple  linear  and  stiff 
systems.  The  Runge  Kutta  routine  was  not  able  to  handle  the 
severe  nonlinearities  when  solving  the  system  of  equations 
in  the  proximity  of  the  digester  heating  zones.  This  was 
apparently  due  to  a  constant  integration  step  size  that 
could  not  be  reduced  in  such  cases.  The  IMSL  routine  which 
has  integration  step  size  control  was  then  chosen  due  to 
this  reason. 

The  steady  state  solution  problem  was  solved  using 
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the  Brown's  modified  Newton-Raphson  method  (IMSL  (39)). 

This  IMSL  routine  was  also  found  to  be  superior  as  it  is  at 
least  guadr at ically  convergent  and  the  number  of  functional 
evaluations  per  iterative  step  are  of  the  order  of  0.5*(N*N 
+  3*N)  as  opposed  to  N*N  for  the  Newton  Raphson  method. 

The  kinetic  model  testing  was  done  by  writing  a 
model  for  a  batch  digester.  This  model  was  basically  the 
same  as  for  the  continuous  case  since  the  same  number  of 
differential  equations  were  to  be  solved.  The  only 
difference  was  the  absence  of  flow  terms  and  adjustments  had 
to  be  made  accordingly.  The  integration  time  was  the  actual 
cooking  time.  For  testing  purposes  actual  laboratory  cook 
data  was  used  and  so  the  heat  balance  equations  could  be 
neglected  as  in  most  of  the  cases  the  data  was  collected 
under  isothermal  conditions  or  the  temperature  profile  was 
explicitly  known.  Data  from  Daleski  (21)  and  from  Smith  and 
William  (88)  was  used.  The  kinetic  parameters  used  for  the 
first  set  of  runs  were  supplied  by  Consolida ted-Bathurst 
Ltd.  (6) ,  but  these  values  had  to  be  tuned  to  get  better 
results.  Seme  of  the  batch  cook  runs  are  shown  in  Appendix 
B.3.  It  was  the  constants  obtained  from  these  runs  that 
were  used  for  subsequent  simulations.  These  constants  are 
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shown  in  Table  2. 

Trivial  testing  of  the  program  was  to  check  the 
overall  working  of  the  program.  Since  individual  routines, 
the  integration  method,  and  the  kinetic  model  were  already 
tested,  it  was  assumed  that  any  other  irregularities  in  the 
model  would  become  apparent  in  this  testing.  The  trivial 
testing  included  such  cases  as  where  the  reaction  rate 
constants  were  set  to  zero  and  the  digester  started  up  from 
a  constant  set  of  initial  conditions.  The  results  from  this 
case  were  that  the  digester  wood  profile  remained  constant 
throughout.  A  similar  test  was  starting  with  a  digester  of 
fully  cooked  chips  and  then  setting  the  reaction  constants 
to  zero.  This  run  showed  the  movement  of  the  cooked  chips 
down  the  digester  and  the  top  steadily  filling  up  with  the 
uncooked  chips.  The  dynamics  of  the  model  were  tested  by  a 
similar  start  up  but  with  the  regular  kinetic  constants  in 
use.  The  dynamics  of  the  model  were  also  tested  by  allowing 
the  digester  to  reach  a  steady  state  value  and  this  value 
compared  with  that  from  the  steady  state  model.  The  results 
were  fairly  close  and,  since  the  same  results  were  obtained 
from  two  independent  programs,  was  a  check  on  the  dynamic 
program.  Since  actual  digester  operating  data  was  lacking  a 
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direct  check  could  not  be  made-  However  it  was  expected 
that  the  model  would  be  able  to  come  close  to  the  actual 
Kappa  number  if  data  were  available  especially  since  cooking 
in  the  wash  zone  was  included- 


CHAPTER  FIVE 


OPEN  LOOP  DIGESTER  SIMULATION  RESULTS 

Several  open  loop  simulation  runs  were  conducted 
using  different  values  for  operating  variables,  such  as 
alkalinity  and  cooking  temperature,  to  show  their  effect  on 
internal  conditiors  in  the  digester  as  well  as  on  the  end 
product.  Most  figures,  e. g.  Figure  5.2,  are  plots  of  the 
state  variable  of  interest  versus  digester  height.  Digester 
height  is  plotted  as  the  vertical  axis  so  that  there  is  a 
direct  correspondence  between  vertical  position  on  the 
figure  and  vertical  position  in  the  actual  digester.  For 
all  the  runs  presented  in  this  chapter,  the  digester  was 
modelled  by  a  total  of  24  sections  and  hence  the  results 
produced  by  the  model  consist  of  24  values  of  each  variable, 
at  each  time  step.  For  a  given  profile,  the  24  points  were 
fit  by  spline  functions  and  plotted  as  a  continuous  curve  so 
that  the  dynamic  changes  in  the  profiles  could  be  followed 
more  easily.  Figures  5.9  and  5.10  are  referred  to  as 
•history  plots*  since  they  show  the  actual  state  of  a  pulp 
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sample  from  the  time  it  entered  the  top  of  the  digester 
until  it  exited  via  the  blow  line.  After  extended  periods 
of  steady  state  operation  the  history  plots  become  equal  to 
the  instantaneous  profiles.  However,  during  transient 
periods  these  plots  show  the  history  of  the  current  exit 
pulp  and  hence  would  be  of  assistance  to  the  operators  of 
downstream  process  units.  The  digester  physical  parameters 
and  operating  conditions  for  run  #1  are  presented  in  full  in 
Appendix  C. 2. 

The  simulation  runs  are  grouped  into  six  sets  and 
are  discussed  in  the  following  sections.  The  sets  are  : 


i) 


ii) 

iii) 

iv) 


v) 

and  vi) 


A  run  from  constant  initial  conditions, 
as  in  the  startup  of  a  digester,  to  a  steady  state, 
alkaline  concentration  changes, 
heater  or  cooking  temperature  changes, 
variation  of  cook  zone  residence  time,  by 
altering  the  production  rate  and  by  varying  the 
length  of  the  cook  zone,  through  changing  quench 
and  extraction  conditions, 
production  and  grade  changes 
a  species  change. 
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bach  set  of  runs  is  discussed  separately  and  the 
overall  results  from  these  runs  are  summarized  in  the  last 
section. 

5. 1  Digester  Run  from  Constant  Initial  Conditions 

to  Steady  State 

This  run  illustrates  the  ability  of  the  model  to 
simulate  large  dynamic  changes.  The  digester  was  assumed  to 
have  constant  initial  profiles  of  fresh  uncooked  chips  and 
the  initial  entrapped  and  free  liquor  concentrations 
throughout  the  digester  were  set  to  those  entering  at  the 
top  of  the  digester.  The  entrapped  and  free  liguor 
temperatures  in  the  different  zones  were  equated  to  the 
temperatures  of  the  external  flows  entering  the  respective 
zones.  A  sample  printout  of  the  initial  profile  is  shown  in 
Appendix  C. 2.  The  digester  reaches  a  steady  state  after 
about  ten  hours  whereas  the  total  residence  time  was 
approximately  six  and  a  half  hours.  This  final  steady  state 
is  the  'base  case'  with  which  other  runs  are  compared  and  is 
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used  as  the  starting  initial  profile  for  the  subsequent 
perturbation  examples- 

In  Figure  5-1,  the  wood  yield  profiles  are  drawn. 
The  first  profile  at  time  zero  is  constant  at  one  hundred 
percent,  since  no  wood  is  assumed  to  have  reacted  at  the 
start.  The  first  or  impregnation  zone  was  at  a  temperature 
of  90  C  and  the  entrapped  liquor  concentrations  are  small 
hence  there  is  little  reaction  in  this  zoce  and  the  steady 
state  is  reached  in  less  than  an  hour-  Most  of  the  yield 
loss  at  this  stage  is  for  the  easily  removable  wood 
components  such  as  acetyls,  etc-  The  heater  and  cook  zones 
have  the  temperature  of  the  upper  and  lower  heater  streams 
and  the  wood  reacts  rapidly  as  most  of  the  extractive  phase 
wood  components  are  removed  within  the  first  five  minutes. 

In  the  extraction  zone  the  cooking  temperature  is  reduced 
and  the  spent  liquor  removed,  which  cuts  off  the  reaction 
substantially  and  this  is  clearly  apparent  from  the  first 
few  profiles-  As  cooking  is  allowed  to  proceed  into  the 
wash  zone,  the  reaction  causes  the  wood  yield  to  decrease 
but  not  as  much  as  in  the  cooking  zones  since  the 
temperatures  and  chemical  concentrations  are  lower.  The 
profiles  at  later  times  show  a  lower  yield  than  before  which 
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:  Wood  Yield  Profiles  for  Run  #1  showing 

25  profiles  from  Startup  to  Final  Steady  State 
(The  profiles  are  at  1,  5,  10,  15,  20,  30,  50, 
75,  100,  150,  200,  250,  300,  350,  400,  450, 

500,  550,  600,  650,  700,  750,  800,  850  and 
900  minutes  respectively) 


Figure  5-  1 
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upper  zones. 


is  due  to  the  movement  of  cooked  chips  from 

The  wood  component  profiles  for  lignin,  mannan  and 
xylan,  figures  5.  2  -  5.4  respectively,  follow  the  same 
general  pattern  as  the  yield  profiles.  The  initial 
concentration  of  100$  was  26.6  kg/100  kg  ODW  (Oven  Dry  Wood) 
for  lignin,  16.2  kg/100  kg  ODW  for  mannan  and  9.5  kg/100  kg 
ODW  for  xylan.  The  residual  percentage  is  higher  for  mannan 
and  xylan  due  to  their  lower  reaction  rates.  The  mannan 
extractive  amount  is  about  55%  and  is  larger  than  the  rest 
so  most  of  the  mannan  reacts  by  the  time  the  chips  pass  the 
heating  zones.  The  xylan  extractive  content  and  rate  is 
lower  and  comparatively  less  of  it  reacts  in  the  same 
period.  The  final  Kappa  number  at  steady  state  was  26.2  and 
the  yield  56  percent.  These  results  are  within  the  range  of 
typical  industrial  operations. 

The  entrapped  and  free  liquor  hydroxide  profiles 
are  plotted  in  Figures  5.5  and  5.6  respectively.  The 
sulphide  plots  are  omitted  since  the  alkalinity  is  directly 
related  to  both  and  and  the  sulphide  curves  would  be 
similar.  The  entrapped  liquor  concentration  is  15$  of  the 
free  liquor  concentration  at  the  top  of  the  digester. 

Further  down  the  digester  the  concentration  difference 
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Figure  5,2  :  Lignin  Profiles  for  Eun  #1 

(  Initial  lignin  Concentration  26-6  kg/100  kg  ODtf) 
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MANNAN,  PERCENT  OF  ORIGINAL 


Figure  5.3  :  Mannan  Profiles  for  Eun  #1 

{  Initial  Mannan  Concentration  16.2  kg/100  kg  ODW  ) 
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Figure  5.4  :  Xylan  Profiles  for  Run  #1 

(  Initial  Xylan  Concentration  9.5  kg/100  kg  0DW  ) 
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between  the  two  phases  gets  smaller  and  the  mass  transfer 
rate  decreases.  The  entrapped  concentration  attains  a 
maximum  value  after  one  hour  and  then  decreases  until  a 
steady  state  value  is  reached.  To  clearly  demonstrate  the 
transients  that  are  involved,  the  profiles  are  split  up  and 
the  ones  from  the  initial  to  the  maximum  are  plotted  on  the 
left  half  of  Figure  5.5  and  the  rest  on  the  right.  The  free 
liquor  profiles  on  the  contrary  have  the  maximum  value  at 
the  start  and  the  concentration  profiles  shew  a  decrease 
until  a  steady  state  is  reached.  The  two  sets  of  liquor 
concentration  profiles  demonstrate  the  combined  effect  of 
reaction,  chemical  usage,  diffusion  mass  transfer  in  the 
chip  and  washing,  along  different  points  on  the  curves. 

Since  in  the  early  stages  of  the  cook,  diffusion  is 
the  controlling  factor  rather  than  reaction  rate,  the 
entrapped  liquor  concentration  reaches  a  maximum  value  at 
about  sixty  minutes  whereas  the  free  liquor  concentration 
falls  to  half  its  original  value  in  the  same  time.  The 
reaction  rate  from  here  onwards  is  controlling  and  the 
entrapped  liquor  concentration  eventually  decreases  to  a 
steady  state.  The  steady  state  entrapped  liquor  curve 
demonstrates  the  predominance  of  diffusion  in  the 
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Figure  5. 5  :  Entrapped  Liquor  Hydroxide  Profiles  for 

Hun  #  1 
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FREE  LIQUOR  HYDROXIDE,  KG/CU  M 
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Figure  5. 6  : 


Free  Liguor  Hydroxide  Profiles  for  Hun  #1 


100 


impregnation  zone  as  the  concentration  reaches  a  maximum 
value  close  to  the  upper  heater  zone.  The  concentration 
below  the  upper  heater  zone  falls  steadily  as  reaction 
becomes  the  controlling  factor.  At  the  guench/extraction 
screens,  there  is  a  further  drop  due  to  mixing  of  the 
cocurrent  spent  and  quench  liquor  and  the  countercurrent 
wash  flow.  A  clean  water  wash  is  used  and  the  entrapped 
chemical  reduces  substantially  due  to  washing,  before  the 
chips  leave  the  digester. 

The  free  liquor  hydroxide  profiles  decrease 
throughout  the  digester  as  the  chemicals  diffuse  into  the 
entrapped  liquor.  The  rate  of  transfer  is  dependent  on  the 
concentration  differential  and  hence  the  reaction  rate.  The 
impregnation  zone  has  little  reaction  consequently  there  is 
only  a  slight  decrease  until  the  upper  heater  zone  where 
there  is  a  sudden  decrease.  The  concentration  decreases 
further  in  the  cook  zone.  The  higher  values  in  the  wash 
zone  at  earlier  times  are  due  to  the  initial  concentration 
profile  in  the  zone.  The  movement  of  countercurrent  wash 
liquor  from  the  bottom  which  displaces  some  of  the  free 
liquor  is  mere  apparent  in  the  earlier  profiles. 

The  total  dissolved  solids  or  reaction  products  in 
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the  entrapped  liquor  are  shown  in  Figure  5.7.  These  curves 
move  in  the  opposite  direction  to  the  chemical  reactant 
curves  and  a  maximum  value  is  reached  at  the 
quench/extraction  zones.  This  value  decreases  further  down 
the  digester  as  reaction  products  are  washed  out  from  the 
chips.  The  effect  of  the  wash  heater  near  the  digester  exit 
is  noticafcle  in  tne  curves  with  a  small  bend  at  that  point. 
This  is  due  to  a  higher  mass  transfer  coefficient  which 
causes  an  improvement  in  the  washing  efficiency. 

The  wood  chip  and  entrapped  liquor  temperature, 
upon  which  the  reaction  rate  is  dependent,  is  plotted  in 
Figure  5.Q.  The  temperature  reaches  steady  state  first  in 
the  impregnation  zone  since  there  is  limited  reaction  and  no 
external  heating  flows  and  the  zone  temperature  is  dependent 
solely  on  the  entering  chip  and  liquor  temperatures.  On 
leaving  the  impregnation  zone  the  chips  rapidly  attain  the 
temperature  of  the  upper  and  lower  heater  streams.  Then  as 
the  chips  pass  through  the  cook  zone  their  temperature  rises 
by  another  six  degrees  from  the  exothermic  heat  of  reaction. 
The  quench  flow  produces  a  rapid  drop  in  this  temperature 
but  this  is  still  higher  than  the  temperature  in  the  wash 
zone.  The  temperature  effect  of  the  cold  blow  and  the  wash 
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Figure  5. 1  : 


Entrapped  Liquor  Solids  Profiles  for  Bun  #1 
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heater  at  the  exit  is  also  apparent  in  the  curves. 

Ihe  history*  curves  for  yield  and  lignin  are 
plotted  in  Figures  5.9  and  5,10,  these  curves  show  the 
cooking  conditions  of  the  pulp  leaving  the  digester  at  900 
minutes  during  its  passage  through  the  digester.  The 
residence  time  for  this  set  was  about  400  minutes  and  the 
yield  and  lignin  curves  exhibit  the  same  behavior  as  the 
earlier  wood  component  profiles.  Figures  5.1  and  5-2.  There 
is  some  reaction  in  the  impregnation  zone  followed  by  a 
steady  drop  through  the  heater  and  cooking  zones.  The 
reaction  rate  decreases  sharply  at  the  quench  zone  with  a 
nearly  zero  rate  in  the  wash  zone,  although  the  slight 
decrease  in  yield  and  lignin  is  still  visible.  The  yield 
history  curve  when  compared  with  the  steady  state  profile  of 
the  same  run  (Figure  5-1)  match  closely  except  in  the 
impregnation  zone  since  the  yield  history  curve  includes 
only  the  major  wood  components  and  excludes  the  fast 
reacting  acetyls,  etc.  while  the  instantaneous  profiles 
calculate  the  ‘true'  yield- 
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figure  5. 8  :  #ood  cnip  Teaperature  Profiles  for  Run  # 
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Figure  5. 9 


Wood  Yield  History  for  Wood  Exiting  at 
900  minutes  for  Bun  #1 
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figure  5. 10  :  Hood  Lignin  History  for  Hood  Exiting  at 

900  minutes  for  Run  #1 
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Alkali  Concentration  Change 

An  alkalinity  change  in  the  cooking  liguor  at  the 
top  of  the  digester  is  generally  directly  dependent  on  the 
operation  of  the  chemical  recovery  cycle.  Although  this 
change  can  also  he  accomplished  by  varying  the  black  liguor 
flow,  this  has  its  limitations  since  the  black  liguor  flow 
is  usually  zero.  The  spent  cooking  liguor  is  recausticized 
to  white  liguor  in  the  chemical  recovery  cycle  and  chemical 
losses  are  usually  made  up  in  the  early  part  of  the  cycle. 

A  step  change  in  the  alkalinity  is  hence  usually  made  by 
direct  hydroxide  addition  from  a  storage  tank. 

Step  changes  of  +  20%  from  the  base  steady  state 
case  were  made  to  the  white  liguor  alkalinity  at  the  top  of 
the  digester.  The  ultimate  response  to  this  variation  is  a 
different  rate  of  reaction  and  this  effect  is  apparent  from 
the  blow  line  yield  and  the  lignin  content. 

Figures  5.11  to  5.14  show  the  results  for  the 
alkali  plus  change.  The  yield  at  the  blow  line  falls  to  54% 
compared  to  56%  for  the  base  case  (Run  #1)  and  the  blow  line 
Kappa  number  to  24.0  as  opposed  to  26-2.  The  steady  state 
alkalinity  in  the  chips  increases  and  so  does  the  residual 
solids  at  the  extraction  screen  and  in  the  final  pulp. 

For  the  negative  step  in  alkali  only  yield  and 
lignin  concentrations  are  plotted  as  shown  in  Figures  5.15 
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and  5-16.  The  response  for  this  case  is  in  the  reverse 
direction  to  that  for  the  alkali  plus  change,  as  is 
expected.  The  yield  and  Kappa  number  values  are  higher  and 
58%  and  29. 1  respectively.  The  yield  results  obtained  from 
the  two  alkali  step  changes  give  an  estimate  of  the 
nonlinearities  involved  and  also  demonstrate  that  the 
H-Factor,  which  is  a  function  of  only  cooking  time  and 
temperature  would  would  not  be  an  adequate  variable  for  the 
control  of  the  Kappa  number  in  such  cases. 

5.3  Cooking  Temperature  Change 

The  main  manipulated  variable  for  the  final  product 
quality  in  the  operation  of  the  Kamyr  digester  in  most  mills 
is  the  cook  zone  temperature  as  set  by  the  lower  heater 
temperature.  Increments  as  small  as  a  quarter  of  a  degree 
Celsius  are  used  to  adjust  the  digester  operation.  To  study 
the  effect  of  temperature,  a  step  change  of  *  5  °C  was  given 
to  both  the  upper  and  lower  heater  temperatures  from  the 
base  case.  As  with  the  alkali  change,  the  heater 
temperature  change  ultimately  increases  the  rate  of  reaction 
and  the  lignin  and  yield  move  in  the  same  direction  as  with 
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figure  5.11  :  Rood  Yield  Profiles  for  Pun  #2 , 

a  +20%  Change  in  Alkali  relative  to  the 
Base  Case  (  Run  #1  ) 

(  The  profiles  are  at  1,  20,  50,  100,  150,  200,  250,  300, 
350,  400,  450  and  500  minutes  respectively  ) 
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Figure  5.  12  :  Lignin  Profiles  for  Sun  #2 

(  *20%  Alkali  Change  ) 
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ENTRAPPED  LIQUOR  HYDROXIDE.  KG/CU  M 


Figure  5. 13  :  Entrapped  Liquor  Hydroxide  Profiles  for 

Run  #2  (  +  20%  Alkali  Change  ) 
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ENTRAPPED  LIQUOR  SOLIDS,  KG/CU  M 
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Figure  5. 14  :  Entrapped  Liquor  Solids  Profiles  for 

(+20%  Alkali  Change  ) 
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Figure  5,15  :  Wood  Yield  Profiles  for  Hun  #3, 

a  -20&  Alkali  Change  Relative  to  the 
Base  Case  {  Pun  #1)  The  profiles  are  plotted 
for  the  same  times  as  in  Run  #2 
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Figure  5.  16  :  Lignin  Profiles  for  Run  #3 

(-20%  Alkali  Change  ) 
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the  alkali  changes.  However,  in  this  case  the  effect 
appears  more  quickly  in  the  digester  as  the  heater  zones  act 
below  the  impregnation  zone  whereas  the  alkali  is  introduced 
at  the  top  of  the  digester. 

The  results  from  the  step  increase  in  temperature 
are  shown  in  Figures  5.17  to  5.20.  The  yield  and  lignin  are 
both  lower  than  the  starting  base  case  as  would  be  expected 
from  a  higher  rate  of  reaction.  The  profiles  for  yield  and 
lignin  (Figures  5.17  and  5.18)  are  similar  to  the  ones  for 
Eun  #1  as  the  heater  zones  are  the  first  affected  and  the 
temperature  change  then  makes  its  way  through  the  lower 
zones.  The  impregnation  zone  conditions  are  unchanged  as 
the  temperature  of  the  streams  entering  the  digester  are  the 
same  as  for  the  base  case.  The  alkali  concentration  (Figure 
5.  19) ,  which  is  initially  the  same  as  the  steady  state  value 
for  Hun  #1,  is  lower  in  the  final  steady  state,  as  more 
chemicals  have  been  used  up  to  dissolve  the  wood  components. 
The  temperature  response  is  plotted  in  Figure  5.20.  For  the 
5  °C  step  down  the  reaction  rates  are  lower  and  hence  more 
wood  and  lignin  appear  in  the  blow  line.  These  results  are 
plotted  in  Figures  5.21  and  5.22  respectively.  The  step 
change  in  the  heater  temperatures  gave  a  Kappa  number  and 
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Figure  5.  17  :  Hood  yield  Profiles  for  Hun  #4, 

a  +5  degree  C  Change  in  Temperature 
Relative  to  the  Base  Case 

(  The  profiles  are  at  1,  5,  10,  15,  20,  25,  30,  40,  50, 
75,  100,  125,  150,  175,  200,  300,  350,  400,  450  and 

500  minutes  respectively  ) 
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figure  5.  18  :  Lignin  Profiles  for  Run  #4 

{  *5  degree  C  Temperature  Change  ) 
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ENTRAPPED  LIQUOR  HYDROXIDE,  KG/CU  M 


Figure  5.  19  :  Entrapped  Liquor  Hydroxide  Profiles  for 

Bun  #4,  (  +5  degree  C  Temperature  Change  ) 
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2  Wood  Chip  Temperature  Profiles  for  Run  #4 
(  +5  degree  C  Temperature  Change  ) 


Figure  5.20 
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wood  yield  of  22.6  and  52%  for  the  plus  change  and  values  of 
31.0  and  60%  for  the  temperature  minus  change.  These 
temperature  changes  are  only  about  +  2%  of  the  cooking  zone 
temperature  but  the  difference  in  the  yield  and  lignin 
content  is  higher  than  that  for  the  alkali  changes  which 
were  approximately  +  20%.  This  provides  sufficient 
justification  for  using  temperature  to  control  the  final 
pulp  quality  in  industry  since  chemicals  are  more  expensive 
and  have  less  effect  on  the  final  product  quality. 

5.4  Cooking  Zone  Length  Variation 

The  variation  of  cooking  zone  length  has  the  same 
effect  as  changing  the  retention  time  in  the  cook  zone. 

This  may  have  to  be  done  in  cases  where  the  pulp  is  not 
fully  cooked  in  the  cook  zone  or  where  the  pulp 
specifications  have  to  be  changed  without  changing  either 
the  production  rate,  temperature  or  the  alkalinity.  The 
length  of  the  cook  zone  would  be  shortest  for  a  ‘hard*  cook 
where  the  upper  heater  is  shut  off  and  the  lower  heater  is 
set  to  the  maximum  possible  under  the  circumstances.  The 
pulp  is  then  given  a  cold  quench  that  cuts  off  the  reaction 
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HOOD  YIELD,  PERCENT 


:  Wood  Yield  Profiles  for  Bun  #5, 

a  -5  degree  C  Temperature  Change  Relative 
to  the  Base  Case  {  Run  #1  ) 

The  profiles  are  at  the  same  times  as  in  Run  #4 


Figure  5.21 
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Figure  5.22  :  Lignin  Profiles  for  Run  #5 

{  -5  degree  C  Teaperature  Change  ) 
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completely  at  the  quench/extraction  zones.  Similarly  the 
length  would  be  longest  for  a  'soft*  cook  where  both  heaters 
are  operating  at  below  their  capacity  and  there  is  little 
quenching,  which  allows  the  reaction  to  proceed  into  the 
wash  zone.  In  this  case  the  length  of  the  cook  zone  is 
dependent  upon  the  extraction  rate  and  the  washing 
temperature . 

The  ’hard1  cook  was  illustrated  by  shutting  off  the 
upper  heater  and  raising  the  temperature  of  the  lower  heater 
by  15  °C.  The  chips  were  allowed  to  cook  at  this 
temperature  and  quenched  with  a  stream  at  about  the  same 
temperature  as  for  the  cold  blow.  This  was  done  by  not 
recycling  quench  liquor  from  the  extraction  screen  and  using 
only  the  'cheater*  flow.  This  cut  off  the  reaction  almost 
completely  and  hence  there  was  limited  cooking  in  the  wash 
zone.  The  yield  and  lignin  curves  are  given  in  figures  5.23 
and  5.24  respectively.  The  yield  and  Kappa  numbers  in  the 
blow  line  were  58%  and  28.0  .  The  entrapped  liquor 
hydroxide  and  the  temperature  are  plotted  in  figures  5.25 
and  5.26. 

The  soft  cook  was  with  the  upper  and  lower  heaters 
at  ten  degrees  below  the  base  level,  and  the  quench  flow 


. 

* 

' 


DIGESTER  HEIGHT,  li 

10.00  20.00  30.00  40.00  50.00  60.00  70. 00 


124 


Figure  5.23  :  Wood  Yield  Profiles  for  Run  #6, 

a  Hard  Cook  of  Chips,  High  Lower  Heater 
Temperature  and  Cold  Quench 

(  These  profiles  are  at  1,  10,  20,  30,  50,  75,  100, 

125,  150,  175,  200,  225  and  250  minutes  respectively  ) 
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Figure  5. 24  :  Lignin  Profiles  for  Run  *6, 

(  Hard  Cook  of  Chips  ) 
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ENTRAPPED  LIQUOR  HYDROXIDE.  KG/CU  M 


Figure  5. 25  :  Entrapped  Liquor  Hydroxide  Profiles  for 

Run  #6,  (  Hard  Cook  of  Chips  ) 
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WOOD  CHIP  TEMPERATURE ,  DEGREES  C 
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Figure  5, 26  :  Chip  Temperature  Profiles  for  Run  #6 

(  Hard  Cook  of  Chips  ) 
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shut  off  completely  to  allow  cooking  to  proceed  into  the 
wash  zone-  The  results  for  yield,  lignin,  entrapped  liquor 
hydroxide  and  chip  temperature  are  shown  in  Figures  5.27  to 
5-30.  The  yield  for  this  case  is  62%  and  the  Kappa  number 
is  32.0  .  The  concentration  of  entrapped  liquor  in  the  pulp 
is  higher  than  normal  and  further  washing  would  have  to  be 
done  before  the  pulp  could  be  used  further  downstream  in  the 
mill. 

The  results  from  the  cook  zone  length  variation 
runs  show  that  the  pulp  quality  is  dependent  on  the  cooking 
temperature  profile  that  the  chips  have  been  subjected  to. 
The  results  also  show  that  although  the  Kappa  numbers  are 
approximately  equal,  the  yields  are  not  and  neither  are  the 
Kappa  numbers.  The  K-number  control  scheme  and  the  H-Factor 
scheme  would  both  have  difficulties  in  predicting  the  final 
product  specifications. 
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Figure  5.27  :  Wood  Yield  Profiles  for  Run  #7, 

a  Soft  Cook  of  Chips,  Low  Heater 
Temperatures  and  Warm  Quench  to  Allow 
Cooking  in  the  Wash  Zone 

These  profiles  are  at  the  same  times  as  in  Run  #6 
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Figure  5.28  z  Lignin  Profiles  for  Fun  11, 

(  Soft  Cook  of  Chips  ) 
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Figure  5.29  :  Entrapped  Liquor  Hydroxide  Profiles  for 

Run  #  7,  {  Soft  Cook  of  Chips  ) 
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figure  5,30  :  Chip  Temperature  Profiles  for  Run  #7 

(  Soft  Cook  of  Chips  ) 
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5.5  Production  and  Grade  Changes 

The  production  rate  and  grade  changes  have  a 
requirement  for  more  or  less  chemicals  and/or  heater 
temperatures  depending  on  the  direction  of  the  change.  The 
Kappa  number  and  yield  would  therefore  be  significantly 
different  if  the  conditions  were  left  unadjusted.  The 
production  change  creates  a  different  cock  zone  retention 
time  and  either  of  the  three  alternatives  discussed  earlier 
could  be  used  to  adjust  the  final  pulp  quality.  A  grade 
change  is  different  than  a  production  rate  change  since  the 
product  specifications  have  changed  and  the  object  is 
minimum  off-grade  pulp.  A  production  rate  change  was  tried 
by  making  a  +10$  step  change  in  chip  flow  rate  relative  to 
the  base  condition.  The  results  for  yield  and  lignin  are 
shown  in  Figures  5.32  and  5.33  and  have,  as  expected,  a 
higher  pulp  yield  and  Kappa  number. 

5.6  Species  Change 


A  species  change  is  accompanied  by  a  different  bulk 
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density  and  different  reaction  rate  constants.  A  species 
change  frcm  hardwood  to  softwood  was  tried  from  base  case 
conditions  and  the  pulp  tracked  through  the  digester.  The 
lignin  and  yield  transient  profiles  are  shown  in  Figures 
5.33  and  5.34.  Hardwood  has  a  faster  reaction  rate  than 
softwood  and  at  steady  state  a  lower  yield  and  Kappa  number 
value  is  observed. 
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Figure  5.31  :  Hood  Yield  Profiles  for  Run  #8, 

a  ♦10X  Change  ic  the  Production  Rate 
Relative  to  the  Base  Case  (  Run  #1  ) 

(  The  profiles  are  at  1,  10,  20,  50,  100,  150,  200, 
250,  300,  350,  400,  450,  500,  550  and  600  minutes  ) 


DIGESTER  HEIGHT.  M 

.00  10.00  20-00  30.00  40-00  50.00  60.00  70.00 


136 


LIGNIN,  PERCENT  OF  ORIGINRL 


Figure  5.32 


Lignin  Profiles  for  Run  #8 
(  *10%  Change  in  the  Production  Rate) 
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WOOD  YIELD,  PERCENT 


Figure  5.33  :  Wood  Yield  Profiles  for  Run  #9, 

a  Species  Change  from  a  Hardwood  to  Softwood 
at  the  Base  Conditions  {  Run  #1) 

The  profiles  are  at  the  same  times  as  for  Run  #8 
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Figure  5.34  :  Lignin  Profiles  for  Run  #9 

(  Species  Change  from  a  Hardwood  to  a  Softwood  ) 
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5-7  Conclusions 

The  main  operating  variables  used  in  digester 
control  were  varied  and  their  effect  on  pulp  yield  and 
lignin  quality  analysed-  As  was  expected,  the  effect  of 
temperature  was  largest  and,  for  a  percentage  change  of  the 
same  magnitude,  the  effect  of  alkalinity  is  the  least-  The 
main  object  of  these  runs  was  to  be  able  to  show  that  the 
scalar  variables  that  are  used  in  the  control  of  product 
quality  were  not  always  unique  and  this  is  demonstrated  by 
the  runs  in  sections  5-2  and  5-4.  These  runs  also  helped  to 
show  that  there  are  a  number  of  alternate  operating 
conditions  that  could  be  used  when  operating  constraints, 
such  as  would  arise  if  a  heater  scaling  problem  was 
encountered.  The  study  of  cook  zone  length  variation  was 
possible  since  the  model  has  the  ability  to  simulate  cooking 
in  the  wash  zone.  The  cook  zone  length  was  shown  to  be  an 
operating  alternative  although  it  is  seldom  used  in  industry 
at  present- 
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CHAPTER  SIX 


CLOSED  LOOP  SCHEMES 

This  chapter  includes  the  review  of  digester 
control  schemes.  First,  digester  control  studies  are 
examined  followed  fay  a  discussion  of  typical  digester 
control  loops  that  are  implemented  in  industry. 

6.1  Review  of  Digester  Control  Schemes 

The  purpose  of  process  control  is  simply  to 
maximize  the  efficient  usage  of  raw  materials  and  process 
equipment.  The  digester  being  at  the  start  of  the  pulping 
process  has  to  be  carefully  controlled  since  production 
upsets  and  disturbances  can  carry  over  and  affect  downstream 
process  units  such  as  the  bleach  plant  and  the  paper 
machines.  Therefore  the  main  objective  in  digester  control 
is  to  maximize  production  by  improving  pulp  yield  and  to 
optimize  the  performance  of  the  digester  so  as  to  minimize 
raw  material,  chemical  and  energy  costs.  The  other  goals 
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are  minimum  variation  in  pulp  quality  as  measured  by  Kappa 
Number,  Permanganate  Number,  Chlorine  Number,  or  fibre 
strength  etc.;  stable  consistent  operation  of  the  process; 
changes  in  production  rate  without  upsetting  pulp  quality; 
changes  of  pulp  grade  or  wood  chip  species  with  minimum 
off-grade  pulp;  and  optimization  of  the  cooking  chemicals 
recovery  cycle  load  by  better  wash  zone  control.  In  recent 
times  with  computer  control  getting  more  popular,  pulp  mills 
are  using  computers  to  give  improved  performance.  The 
reason  for  this  is  more  consistent  plant  operation 
regardless  of  operating  personnel  capabilities. 

Process  control  schemes  and  strategies  in  the  pulp 
and  paper  industry  have  been  reviewed  by  by  Holm  and  Perry 
(42),  Williams  (104)  and  more  recently  by  Perron  and  Ramaz 
(80).  These  articles  report  on  the  philosophy  behind  the 
design  and  implementation  of  control  systems  in  the  industry 
and  have  references  to  mathematical  models  of  the  processes. 

Computer  control  was  introduced  to  the  pulp  and 
paper  industry  in  1961  and  the  first  continuous  digester  was 
brought  under  computer  control  in  1964.  Both  Fox  (28)  and 
Smith  (90)  report  on  the  performance  of  the  first  closed 
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loop  process  computer  control  system  for  a  digester.  There 
are  survey  articles  on  the  extent  of  digital  computer  usage 
in  the  industry  by  Brewster  and  Bjerring  (11),  Keyes  (54) 
and  Gee  amd  Chamberlain  (32)  .  The  Brewster  article  was 
published  when  computer  control  was  still  in  its  inception 
in  the  industry  and  as  of  March  1969,  there  were  a  total  of 
47  computer  installations  in  the  U.  S.  Pulp  and  Paper 
Industry.  In  1976  there  were  over  600  installations  (54) 
including  approximately  33  continuous  digesters  under 
computer  control.  This  rise  in  the  number  of  computer 
installations  is  due  to  the  development  of  control  packages 
by  computer  vendors,  improvement  in  computer  performance, 
decrease  in  hardware  cost  and  the  development  of  sensors 
which  make  computer  control  more  convenient.  Economic 
benefit  from  installing  computers  have  been  reported  (10, 

43,  98)  and  payout  periods  for  such  systems  range  between 
one  and  two  years. 

The  problems  involved  in  measurement,  modelling  and 
control  have  been  described  by  Mayer  (7  3)  and  Church  (20). 
Measurement  problems  are  attributed  to  unreliable  and 
inaccurate  sensors  for  wood  variables  like  chip  moisture 
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content,  wood  quality,  lignin  concentration  and  on  line 
alkali  measurement  for  cooking  chemicals,  A  paper  by 
Macleod  (71)  mentions  that  vendors  are  developing  such 
instruments  and  that  they  are  in  the  course  of  being 
throughly  tested.  Long  time  delays,  measurement  noise  and 
inaccurate  models  in  addition  to  sensors  are  blamed  for  most 
of  the  control  problems. 

Pulping  models  have  been  developed  to  accurately 
describe  the  cooking  process  and  they  have  been  discusssed 
in  an  earlier  chapter.  The  object  of  such  models  has  always 
been  for  an  improved  understanding  of  the  problem  with  an 
ultimate  use  for  better  control  purposes.  The  aim  is  to  be 
able  to  accurately  predict  the  output  controlled  variables 
from  the  process  input  variables.  The  increasing 
utilization  of  computers  for  control  have  made  it  imperative 
for  better  models  so  that  a  feedforward  compensation  type 
decision  can  be  made  before  the  disturbance  significantly 
affects  the  process.  One  of  the  first  studies  in  this 
regard  was  by  Vroom  (99)  who  derived  the  H-factor.  This 
variable  is  used  as  a  combined  measure  of  the  cooking 
temperature  and  the  duration  of  cooking  and  either  of  the 
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two  can  be  adjusted  for  a  change  in  the  other  to  give  a 
consistent  degree  of  cooking.  The  work  followed  by  Carroll 
(15),  Kerr  (52,  53),  McKibbins  et  al  (75),  Service  (86), 
Hatton  (36,  37),  and  Boyle  and  Tobias  (9)  can  be  classed  as 
statistical  or  mathematical  models.  Carroll  (15)  and  Boyle 
and  Tobias  (9)  have  models  for  the  overall  pulping  process 
while  the  others  are  limited  to  digester  operation.  These 
models  relate  the  output  variables  like  pulp  yield  or  Kappa 
number  to  input  variables  such  as  cooking  time,  temperature, 
alkali  to  wood  ratio  and  H-  Factor  for  a  particular  wood 
species.  The  more  detailed  models  include  chip  moisture  and 
size  and  cooking  liquor  sulphidity.  Other  more  detailed 
models  by  Stillman  (93),  Johnsson(49)  and  Smith  and  Williams 
(88)  have  been  described  earlier. 

The  increasing  use  of  computer  control  has  involved 
computer  vendors  in  the  development  of  various  models  and 
control  packages  for  pulp  and  paper  applications  including 
continuous  digester  control.  These  packages  are  available 
on  a  turnkey  basis  and  Church  (20)  and  Gee  and  Chamberlain 
(32)  have  lists  of  vendors  and  their  applications  by  process 
area.  The  Kappa  Number  control  loop  is  considered  the  most 
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complex  as  it  employs  a  model  of  the  kinetics.  Measurex 
(10,  64)  use  their  batch  digester  model  and  a  model  of  the 
pulp  movement  from  the  heating  zone  to  the  blow  valve. 
Industrial  Nucleonics  (30)  use  an  H-Factor  model 
incorporating  the  lower  cooking  zone  volume.  This  model 
assumes  that  the  initial  alkali  concentration  is  constant 
(20).  Foxboro  (4)  uses  a  model  reference  control  technique 
with  a  Kappa  number  steady  state  model  based  on  the 
Alkali/Wood  ratio,  chip  temperature  and  the  residence  time 
and  a  dynamic  model  which  includes  dead  time  compensation  on 
the  output  from  the  static  model  before  the  resultant 
predicted  valve  is  compared  with  a  laboratory  measured  Kappa 
number  value.  The  error  is  used  to  adjust  parameter  valves 
in  the  static  and  dynamic  models.  Now  some  organizations 
are  buying  their  own  hardware  and  developing  software  "in 
house"  to  match  their  particular  need.  "Modern"  control 
techniques  have  been  mostly  used  on  paper  machines  in  the 
past  but  not  for  continuous  digesters.  Cegrell  and  Hedguist 
(18)  have  used  linear  stochastic  control  theory  to  derive  a 
self-tuning  regulator  algorithm  to  control  the  Kappa  number. 

A  survey  conducted  on  pulping  methods  of  125  mills 


!  . , 


' 

♦  I-  \i  s  >  t  i  0  5  1  1  ’*■  ;  f  "••  •’  *  I 


:  ,  -,  •  M  u.ii  tac-:  u 


146 


in  North  America  (3Sf  39)  reported  that  905?  of  the  mills 
considered  the  cooking  temperature  as  a  valuable  control 
variable  to  correct  the  control  of  a  cook,  44%  of  them  used 
alkali  concentration  and  27%  used  residence  time  in  the  cook 
zone  while  25%  used  all  three,  A  digester  model  in  a  three 
variable  situation  would  therefore  make  possible  an 
economically  optimum  adjustment  among  the  variables 
considering  pulp  production  and  mill  restraints  on  chemicals 
and  energy. 

The  digester  control  system  and  loops  are  shown  in 
the  piping  and  instrumentation  diagram.  Figure  3, 1  .  The 
control  system  is  hierarchical  in  design-  Level  zero,  the 
lowest  level,  is  the  measurement,  data  acquisition, 
scanning,  alarming  and  the  process  operators  console.  The 
next  higher  level  one  is  single  loop  setpoint  regulatory 
control  for  each  of  the  control  variables.  Level  two  is 
multivariable  direct  control  which  could  include  cascade, 
ratio  or  feedforward.  The  levels  one  and  twc  control  loops 
receive  their  setpoints  from  the  next  higher  level.  Level 
three  is  supervisory  control  for  the  overall  process,  for 
production,  yield,  pulp  quality  and  grade  changes.  The 
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regulatory  control  objectives  for  the  three  levels  are  chip 
feed  rate  to  the  top  of  the  digester,  blow  flow,  and  chip 
level,  chemical  and  liquor  to  wood  ratios,  quench  and 
extraction  temperatures,  wash  zone  countercurrent  flow 
control.  Kappa  Number  and  Yield.  These  loops  are  described 
in  more  detail  in  the  following  sections. 

6.2  Ch ipme te r  Control 

Chipmeter  rate  control.  Figure  6.1  regulates  the 
chipmeter  speed  to  maintain  a  constant  dry  wood  rate  to  the 
digester.  The  chipmeter  is  dynamically  calibrated  each  time 
the  chip  bin  is  emptied,  as  it  is  a  volumetric  device  and 
the  filling  efficiency  and  consequently  the  chipmeter 
capacity  per  revolution  are  required.  The  chip  moisture  for 
an  oven  dry  rate  calculation  is  either  calculated  from  a 
steaming  vessel  heat  and  material  balance  (27)  which  is  not 
very  accurate  and  large  errors  are  not  uncommon  or  it  is 
read  directly  from  an  on-line  sensor.  These  chip  moisture 
meters  will  measure  moisture  content  during  summer  months  in 
Canada  but  are  severely  limited  when  the  chips  are  frozen. 
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Chip  Feed  Rate  and  Level  Control 
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The  chip  packing  density,  moisture  content  and  desired  yield 
being  known  the  computer  is  fed  a  specific  production  rate 
value  and  the  desired  chipmeter  rpm  is  evaluated. 

6.3  Chemical  and  Liquor  to  Wood  Ratios  Control 

The  flow  of  cooking  liquor  to  the  top  of  the 
digester  is  controlled  to  maintain  the  ligucr  to  wood  ratio 
and  the  alkali  to  wood  ratio-  The  liquor  requirements  vary 
between  mills  and  with  the  product  produced-  Typically, 
bleached  papers  require  more  liquor  and  chemicals-  The 
liguor  ratio  could  vary  from  3  to  6  cubic  metres  of 
liquor/ton  of  dry  wood  and  the  amount  of  active  alkali  from 
10  to  20  percent  (kg  of  AA  as  Na20  /  100  kg  of  dry  wood) . 

The  chemical  to  wood  ratio  is  set  to  keep  the 
cooking  chemicals  slightly  in  excess  of  that  required  for 
delignif ication.  The  actual  chemical  input  is  calculated 
from  white  liquor  concentration  and  flow  and  dry  wood  flow 
rate  from  the  chipmeter.  The  white  liquor  flow  is  used  as  a 
manipulated  variable  to  correct  any  imbalances.  This  is  a 
sluggish  control  loop  with  a  time  lag  of  about  one  hour  as 
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Figure  6.2  :  Alkali/Wood  and  Liquor/Wood  Ratio  Control 
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the  concentration  in  the  top  circulation  rises  slowly-  For 
a  rapid  change  the  chemical  concentra tion  is  allowed  to 
overshoot  the  required  value. 

The  liquor  to  wood  ratio  control  is  maintained  by 
manipulating  the  flow  of  weak  black  liquor  to  the  top  of  the 
digester-  The  liquor  downflow  in  the  impregnation  and 
cooking  zones  is  calculated  from  a  volumetric  flow  balance 
about  the  digester-  The  liquor  downflow  originates  from 
chip  moisture,  white  liquor,  steam  condensate  from  the 
steaming  vessel  and  the  top  dilution  black  liquor-  The 
amount  of  dilution  liquor  added  must  be  balanced  against  the 
heat  load  on  the  upper  and  lower  heaters  and  chemical 
recovery  load  at  the  extraction  screens.  The  alkali 
concentration  is  usually  manually  tested  every  hour  and  a 
reliable  outline  alkali  sensor  would  be  an  advantage  for 
both  loops. 

6-4  Chip  Level  Control 

The  chip  level  and  blow  flow  are  interrelated  as 
they  affect  the  chip  residence  time  in  the  digester  zones. 
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The  chip  level  is  evaluated  from  the  electric  current  load 
of  the  stirrer  in  the  digester  top  separator.  Figure  6.1. 

The  chip  level  controller  then  either  adjusts  the  chipmeter 
speed  and  hence  the  chip  flow  at  the  top  of  the  digester  or 
the  blow  flow  rate  and  consistency. 

6.5  Blow  Flow  and  Consistency  Control 

The  blow  flow  and  chip  level  are  controlled  to  give 
a  constant  residence  time  in  the  various  digester  zones  and 
primarily  the  cook  zone  (Figure  6.3).  The  rate  of  fibre 
discharge  is  dependent  on  the  blow  flow  and  its  consistency 
and  are  set  to  the  desired  production  rate.  The  consistency 
is  controlled  by  the  bottom  dilution  flow,  digester  pressure 
at  the  bottom  or  the  bottom  scraper  rpm.  The  fibre 
discharge  rate  and  the  chip  feed  rate  are  in  turn  used  in 
the  control  of  chip  level  in  the  digester. 

6.6  Temperature  Control 


Temperature  of  the  chip  mass  is  considered  the  most 
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Blow  Flow  and  Consistency  Control 
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important  pulping  variable  by  most  mills  (38)  ,  as  was 
mentioned  earlier.  This  is  due  to  the  large  sensitivity  of 
temperature  in  determining  the  final  product  quality  (Kappa 
Number)  of  the  pulp. 

There  are  normally  two  separate  heating  zones  in  a 
Kamyr  digester,  refered  to  as  the  upper  and  lower  heating 
zones.  Only  the  controls  for  the  lower  zone  are  shown  in 
Figure  6.4  as  the  upper  zone  control  is  similar. 

The  heater  flows  and  heater  entrance  and  exit 
temperatures  are  the  only  measurable  variables  for 
temperature  control  of  chips.  Since  a  direct  temperature 
measurement  of  the  pulp  mass  is  not  available,  the  value  is 
estimated  using  a  model  of  the  heater  zones.  This  model  is 
just  a  dynamic  energy  balance  based  on  the  heater  section. 
The  setpoints  for  the  two  heaters  are  set  by  the  Kappa 
number  control  scheme.  The  upper  heater  chip  exit 
temperature  is  estimated,  compared  to  the  setpoint,  and  the 
error  used  as  a  basis  for  feedforward  compensation  in  the 
lower  heater  controller.  The  time  lag  for  a  change  in  the 
two  heater  temperatures  varies  from  five  to  ten  minutes 
depending  mainly  on  the  heater  flow  rate  entering  the 
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digester. 


6.7  Quench  and  Extraction  Control 

After  cooking  is  complete  the  spent  liquor  is 
displaced  and  removed  at  the  screens  before  it  is  sent  for 
recovery.  The  cooking  reaction  is  reduced  by  the  quench 
liquor  just  above  the  extraction  screens  as  shown  in  figure 
6.5.  The  boundary  between  the  effective  cook  zone  and  the 
wash  zone  is  controlled  by  the  extraction  flow  rate  and  the 
temperature  of  the  cooler  countercurrent  flow  of  liquor  from 
the  wash  zone.  A  cheater  line  from  cold  blow  to  the  quench 
line  gives  independent  control  of  the  quench  temperature. 

If  the  cheater  line  is  not  available  then  the  quench 
temperature  is  dependent  only  on  the  extraction  rate. 

The  extraction  flow  determines  the  wash  liquor 
upflow  and  is  manipulated  to  keep  the  wash  flow  at  its  set 
point.  The  countercurrent  wash  flow  is  estimated  in  one  of 
two  ways  ways,  a  heat  balance  on  the  wash  heater  can  be  used 
since  the  heat  load  is  proportional  to  the  upflow  or  through 
a  volumetric  balance  at  the  bottom  of  the  digester.  Wash 
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flow  is  characterised  by  the  dilution  factor  which  is  the 
ratio  of  wash  liquor  upflow  to  the  pulp  downflow  in  the 
washing  zone  and  typical  values  are  betweeen  2  and  3  kg  of 
water/kg  of  air  dry  pulp,  A  smaller  dilution  factor  would 
result  in  insufficient  washing,  i.e.  reduced  washing 
efficiency,  and  a  larger  value  would  put  an  additional 
burden  on  the  wash  heater  and  the  recovery  cycle-  A  large 
dilution  factor  would  also  affect  the  cook  zone  boundary  at 
the  extraction  screens  while  a  small  value  would  extend  the 
cook  zone  almost  to  the  bottom  of  the  digester.  This  is 
sometimes  used  by  the  operators  if  there  is  insufficient 
cooking  in  the  cook  zone. 

6.8  Ka ppa  Number  and  Pulp  Yield  Control 

The  Kappa  Number  and  pulp  yield  are  interrelated 
and  equations  have  been  developed  that  relate  the  two 
variables  for  a  particular  wood  species,  e.g.  Johnsson 
(49),  Eydholm  (81  p635)  and  Smith  and  Williams  (88). 
Therefore,  in  most  cases  the  pulp  yield  versus  the  lignin 
concentration  relationship  is  known  by  mills  for  the  wood 
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species  that  they  process.  The  Kappa  number,  and  hence  the 
yield,  is  usually  regulated  by  the  upper  and  lower  heater 
outlet  temperatures.  Since  the  Kappa  number  is  the  most 
important  variable  for  determining  final  pulp  quality 
precise  temperature  control  is  vital.  The  H-Factor  which 
combines  both  the  cooking  time  and  temperature  is  also  used 
as  a  measure  of  the  degree  of  cooking.  The  aim  is  to  cook 
the  wood  chips  to  a  lignin  content  in  its  bulk  reaction  rate 
phase  instead  of  the  residual  reaction  phase  where  the  rate 
of  cellulose  removal  is  faster  than  that  of  lignin.  The 
temperature  and  the  residence  time  in  the  cook  zone  are  the 
only  two  variables  that  would  affect  the  pulp  quality  once 
the  chips  are  in  the  cook  zone.  The  cook  zone  residence 
time  is  only  adjusted  when  a  change  in  production  rate  is 
desired,  if  there  is  a  restraint  on  steam  flow  and  a  higher 
cook  temperature  is  desired,  or  when  the  chips  are  not 
sufficiently  cooked  and  the  operator  would  like  the  cooking 
to  proceed  into  the  wash  zone. 

A  steady  state  model  of  the  digester  is  used  to 
calculate  the  Kappa  Number  based  on  the  Alkali/Wood  ratio, 
the  cooking  temperature  and  the  residence  time  in  the  cook 
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zone  as  shown  in  Figure  6.6.  One  of  the  three  operating 
variables  is  adjusted  based  on  the  current  operating 
constraints.  The  pulp  yield  or  the  Kappa  Number  from  the 
blow  line  is  analysed  typically  every  hour  and  the  error  fed 
back  to  correct  for  drift  in  the  model  parameters.  The 
model  is  used  to  compensate  for  the  long  time  delay  between 
the  end  of  the  cook  zone  and  the  time,  a  sample  is  analysed 
before  the  lower  heater  temperature  is  adjusted. 

6.9  Production  and  Grade  Change  Control 

Production  rate  and  grade  changes  can  be  made  in  a 
scheduled  or  unscheduled  manner  and  both  have  different 
operating  procedures.  The  planned  changes  can  be  due  to 
production  demand,  a  different  grade  requirement  or 
introduction  cf  a  different  wood  species.  An  unscheduled 
change  is  usually  due  to  a  bottleneck  in  the  chipyard  or 
chip  bin  or  to  meet  sudden  upstream  or  downstream  demands 
for  a  change  in  production  rate.  The  goal  in  both  cases  is 
to  make  the  switch  as  smooth  as  possible  with  minimum 
disturbance  towards  quality  and  yield. 
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For  scheduled  changes  the  heater  temperature  is 
gradually  changed  and  the  chemicals  suddenly  changed  to  a 
new  value  that  is  followed  by  a  gradual  change  in  the 
production  rate.  For  an  emergency  change  the  chip  feed 
rate,  heating  zone  temperatures  and  liquor  flows  are  rapidly 
changed  to  new  setpoint  values.  This  is  to  minimize  the 
pulp  quality  and  yield  transients  which  cannot  otherwise  be 
avoided.  A  grade  change  is  similar  to  a  production  change 
and  is  treated  similarly.  This  change  can  be  due  to  a  yield 
or  Kappa  number  change  with  the  same  species.  A  change  of 
wood  species  results  in  changes  in  bulk  density,  cooking 
rate  and  yield.  The  object  is  to  produce  pulp  which  has  the 
desired  properties  when  the  new  species  arrives  at  the  blow 
line  without  producing  a  hard  or  soft  pulp  during  the 
interim  period.  This  is  achieved  by  changing  the  operating 
parameters  as  the  species  interface  moves  through  various 
digester  zones.  The  parameters  that  are  changed  are  the 
liquor  to  wood  and  chemical  to  wood  ratios,  cooking 
temperature  and  the  wash  target.  An  unscheduled  species 
change  requires  rapid  computation  of  new  conditions  but  they 
are  implemented  in  an  identical  manner. 


' 

•:  U  "  r..  •  -  *  '  -  ' 

■ 


CHAPTER  SEVEN 


CONCLUSIONS  AND  FUTURE  WORK 


The  conclusions  and  recommendations  for  further 
work  are  presented. 

7.1  Conclusions 


A  nonlinear,  dynamic  model  of  a  continuous 
digester  has  been  developed  and  used  in  open  loop  s 
studies.  Since  industrial  digesters  do  not  normail 
extensive  number  of  sampling  ports  or  sensors,  it  h 
been  possible  to  obtain  the  data  necessary  for  a  di 
comparison  of  concentration  profiles  and/or  time  hi 
However,  individual  modules  of  the  program  have  bee 
separately  and  the  total  program  produces  results  t 
intuitively  correct  and  which  compare  favourably  wi 
results  of  previous  studies. 

Some  compromises  and  approximations  were  m 
especially  in  the  area  of  chemical  kinetics  to  keep 
overall  program  size  and  execution  time  to  a  reason 
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level  but  the  program  is  very  general  in  terms  of  what 
configurations,  operating  conditions  and  disturbances  it 
will  handle.  All  the  digester  dimensions,  physical 
properties,  physical  properties  etc.  can  be  set  by 
modifying  input  data,  or  in  some  cases  portions  of  specific 
subroutines.  Operating  conditions  can  be  set  by  the  user 
and  are  specifically  programmed  to  facilitate  closed  loop 
control  studies.  The  program  will  simulate  the  dynamic 
response  to  chip  moisture  disturbances,  species  changes, 
changes  in  inlet  liquor  flow  and/or  concentration,  plus 
changes  in  heater  circulation  flows  and/or  temperatures. 

The  model  includes  thirteen  state  variables  in  the  dynamic 
equations,  allows  for  all  cooking  reactions  to  proceed  into 
the  wash  zone,  and  includes  optional  features  such  as  the 
use  of  ‘cheater*  flow  for  the  independent  adjustment  of 
quench  temperature. 

Separate  versions  of  the  program  are  also  available 
to  simulate  batch  digesters  and  to  produce  steady  state 
profiles  and  histories.  These  programs  were  used  to  produce 
an  'independent*  check  on  results. 

The  results  of  the  general  simulation  program  are 
available  in  three  forms.  First,  as  profiles  of  any  of  the 
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state  or  derived  variables  such  as  lignin  concentrations, 
temperature.  Kappa  number  etc.  plotted  versus  digester 
height.  These  plots  define  the  operation  of  the  digester  at 
any  instant  of  time-  Second,  the  program  can  produce 
'history  plots'  of  operating  conditions  or  variables,  such 
as  lignin  concentration  versus  time.  These  plots  define  the 
operating  conditions  used  and  the  state  variables  which 
characterizes  the  pulp  from  the  time  a  particular  unit  of 
pulp  enters  the  digester  until  it  leaves  the  blow  line. 
Third,  the  program  produces  listing  of  the  steady  state 
profiles  plus  summaries  of  the  operating  conditions  and 
configuration.  One  of  the  recommendations  arising  from  this 
work  is  that  industry  consider  generating  the  profiles 
on-line  and  use  them  for  operator  guidance  or  process 
control.  This  should  permit  significant  improvements  in 
product  quality  and  uniformity  because  at  the  present  time 
the  exit  pulp  quality  cannot  be  measured  directly  by  any 
on-line  set  of  measurements  nor  can  the  inlet  digester 
conditions  required  to  produce  a  pulp  of  specified  quality 
be  easily  determined-  Presently,  the  operators  use  scalar 
measures  of  pulp  quality  such  as  H-Factor,  residual  lignin 
concentration.  Kappa  number  etc.  for  control  purposes. 
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However,  as  shown  by  simulation  results  from  this  study  ( 
and  as  is  undoubtly  found  in  industrial  practice  )  the  same 
scalar  measures,  such  as  Kappa  number  at  the  blow  line,  can 
be  produced  using  significantly  different  operating 
conditions.  It  is  almost  certain  that  the  actual  production 
history  of  a  particular  unit  of  pulp  would  influence  its 
downstream  processing  characteristics  such  as 
1 bleachabil ity '  and/or  final  product  specifications  such  as 
paper  strength,  smoothness  etc.  The  display  of 
instantaneous  digester  profiles  should  make  it  possible  for 
the  digester  operator  to  control  the  digester  operation  more 
closely  and  to  improve  reproducibility.  Use  of  displays 
showing  the  production  history  of  the  unit  of  pulp  currently 
entering  a  downstream  process,  such  as  the  bleach  plant, 
should  improve  control  and  product  quality  in  these 
downstream  units. 

Use  of  a  simulation  program  such  as  the  one 
developed  in  this  work,  in  parallel  with  actual  digester 
operation  would  permit  the  process  operators  to  evaluate 
alternatives  for  implementing  chip  species  changes,  grade 
changes  and/or  production  rate  changes.  Displays  of 
instantaneous  profiles  and  historical  plots  would  show  the 
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progress  of  changes  through  the  digester  and  hence  provide  a 
basis  for  operator  control  to  improve  response  times  and 
reduce  off  grade  pulp.  Careful  tuning  of  the  parameters  of 
the  dynamic  model  would  be  required  but  simulation  of 
startup  dynamics  (c. f .  Figures  5. 1-5.5)  and  optimization 
studies  to  reduce  energy  consumption  or  balance  tradeoffs 
between  chemical  consumption  and  pulp  quality  are  also 
possible  (c.f.  Chapter  5,  where  the  effects  of  changes  in 
cooking  temperature  and  chemical  concentration  were 
demonstrated).  The  parameters  in  the  dynamic  model  such  as 
reaction  rate  constants,  diffusion  coef f icients ,  heat 
transfer  coefficients,  pulp  density  etc.  plus  their 
dependence  on  operating  conditions  such  as  temperature, 
should  obviously  be  ’tuned1  via  a  continuing  project 
involving  production  and  development  personnel. 

7.2  future  Work 

As  suggested  above  the  dynamic  model  developed  as  a 
part  of  this  project  has  significant  potential  for  effective 
industrial  utilisation.  However,  one  of  the  prime 
motivations  in  developing  this  model  was  to  provide  a 
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convenient  and  efficient  means  of  evaluating  improved 
computer  control  techniques  -  especially  at  the  higher  or 
supervisory  levels.  Some  specific  approaches  that  appear  to 
have  significant  potential  for  industrial  use  include  : 

i)  Material  and  energy  balance  control  schemes  to 
coordinate  the  supervision  of  individual  DDC  loops. 

ii)  On-line  constraint  identification,  using 
techniques  such  as  linear  programming  followed  by 
appropriate  modification  of  operating  conditions  and/or 
ob  ject ives. 

iii)  Use  of  modern  multivariable  control  techniques 

such  as  those  investigated  previously  in  this 
department  (26) .  1 

For  example  : 

a)  time  delay  compensation  (38-A) , 

b)  frequency  domain  designs  to  'decouple1,  i.e. 
reduce  interactions  between  controlled  variables 
(62)  , 

1  the  following  references  are  to  the  individual  theses. 
Several  papers  based  on  these  projects  were  published 
subsequently  in  different  journals. 


' 


s 


169 


c)  use  of  filters  and/or  observers  to  estimate 
variables  that  cannot  be  measured  directly  (33, 
79), 

d)  control  of  closed  loop  system  response  and/or 
parameter  sensitivity  by  eigenvalue  and 
eigenvector  assignment  (87) , 

e)  ‘optimal*  manipulation  of  input  variables  to 
achieve  minimal  time  response  in  the  presence  of 
operating  constraints  (77). 

iv)  Self  tuning  regulators  (19,  82),  and  hyperstable 
adaptive  control  systems  that  require  only  input/output 
measurements  for  the  plant  (50). 

v)  Supervisory  and  distributed  system  control 
concepts  that  would  permit  feedforward,  feedback  and 
adaption  of  information  among  all  the  major  units  in  a 
mill  from  the  digester  to  the  paper  machine.  For 
example,  factors  such  as  the  bleac hability  of  the  pulp 
and  the  quality  of  paper  produced  from  it,  provide  the 
ultimate  means  of  evaluating  digester  operations  and 
globally  optimize  economic  behaviour.  A  computer 
system  is  required  for  such  feedback  and  evaluation. 
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vi)  The  kinetic  model  used  in  this  study  is  too 
simple  to  accurately  represent  the  complex  set  of 
pulping  reactions.  Present  research  includes 
development  of  processes  to  reduce  hemicellulose 
degradation  by  addition  of  chemicals  such  as 
anthraguinone  in  the  cooking  liquor.  Studies  have  also 
been  done  for  reprecipitation  of  the  carbohydrates 
during  the  later  portion  of  the  cook.  Accurate 
simulation  models  for  such  processes  are  required. 
Cooking  data  for  various  wood  species  which  would 
assist  in  parameter  estimation  for  these  models  is  also 
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APPENDIX  A 


DERIVATION  OF  MODEL  EQUATIONS 

This  appendix  is  in  support  of  material  covered  in 
chapter  three.  The  partial  differential  equation  is  derived 
first  and  the  set  of  thirteen  nonlinear  ordinary  differential 
equations  in  the  next  section. 


A.l  Partial  Differential  Equation 

The  unsteady  state  wood  component  mass  balance  was 
derived  in  cylindrical  coordinates.  The  equation  can  also  be 
found  in  a  standard  textbook  (eg.  Bird,  Stewart  and  Lightfoot, 
6,  p559,  Table  18.2-2  Eqn  B)  : 


3C 

at 


r  i 

r  2  2  ) 

Jv  3C  vG  SC  v  3C 

13  r  3C  ,  1  a  c  ,  a  c 

r  +6  +  z 

=  D; 

2-  2  „  2 1 

^  3r  r  30  8zJ 

r  3r  3r  r  36  3z  J 

+  R 


(A. 1-1) 


where  C  is  a  wood  component  concentration,  vr»  vQ,  and 

V  are  the  velocity  components  in  the  three  directions,  r  is  the 
z 

radius,  D  is  the  diffusion  coefficient  and  R  the  total  reaction 
term  i.e.  the  rate  of  production  per  unit  volume..  (  The  nomenclature 
is  summarized  in  a  table  at  the  end  of  this  appendix.  ) 
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Neglecting  the  radial  (r) ,  angular  (0 Components  and  the 
diffusion  coefficient,  D.  Equation  (A. 1-1)  reduces  to  : 


8C  ,  v  3C 
+  z  _ 

at  3z 


R 


0 


(A. 1-2) 


The  same  equation  can  be  written  in  functional  terms  as: 


3C(z,t)  +  vz(z)  3C(z,t)  _  R(c>z>t)  .  0  (A. 1-3) 

3 1  3z 


The  boundary  conditions  for  this  case  are  : 


=  C  (t)  at  z=0  and  t>0  (A. 1-4) 

o  = 

=  C(z)  at  t=0  and  CKz^H 


C(0,t) 
C(z ,0) 


(A. 1-5) 
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Discretizing  equation  (A. 1-2)  around  z 


dC  (t) 
n 

dt 


+ 


v  AC  (t) 
z,n  n 


Az 


R  (C  ,t)  =  0 

n  n 


(A. 1-6) 


or  when  multiplying  throughout  by  the  volume  of  section  n  in  the 
digester,  V 

n 


VndCn(t)  +  VnVz,n  Cn-l(t)  -  Cn(t)  -  V  R  (C  .,)  ■  0 

-  n  n  n 


dt 


Az 


or,  assuming  constant  crossectional  area  in  each  section 


Vn  dCn(t)  +  Q  (t)  { C  ,(t)  -  C  (t)>  -  V  R  (C  ,t)  =  0  (A. 1-7) 

-  n  n— i  n  n  n  n 

dt 


The  boundary  conditions  (A. 1-4)  and  (A. 1.5)  reduce  to  : 


C(0,t)  «  C  (t)  at  z=0 

o 


(A. 1-8) 


and 


C(z  ,0)  =  C  at  t=0 

n  n,o 


(A. 1-9) 


respectively,  where  Cq  is  the  concentration  value  entering  at  the 

top  of  the  digester  and  q  is  the  initial,  starting  concentration 

in  section  n,  and  Q  is  volumetric  flow  rate,  i.e.,  Q=V  v  /Az. 

n  n  z,n 
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The  same  equations  if  derived  for  the  chemical  components 
in  the  free  and  entrapped  cooking  liquor  would  have  an  additional 
diffusion  term  for  mass  transfer  between  these  two  phases. 


A. 2  Ordinary  Differential  Equations 

The  set  of  nonlinear  ordinary  differential  equations  are 
presented  here  for  a  digester  section  n  : 

Overall  mass  balance  : 


where 


Vn  dPn  =  V  Func  (z  2  ,  HFAC »  t) 
-  n  n-1  n 


dt 


(A. 2-1) 


V  is  volume  of  section,  m 
n 

3 

is  the  density  (compaction) ,  kg/m 

z  - ,  z  is  position  of  section  n,  m 
n— l  n 

HFAC  is  the  H-Factor  and  a  measure  of  degree  of  cooking 

Individual  wood  component  balances  : 

The  individual  wood  component  balance  in  terms  of 

X  ,  the  wood  specie  concentration  in  section  n  is  : 
n 


p  V  dX 
n  n  _ n 

dt 


w 


n-1 


w  X  -  R 
n  n  x  ,n 


(A. 2-2) 


where 


w  is  the  wood  flow  rate  from  n  in  kg  ODW/min 
n 

and  R  is  the  total  reaction  term,  hydroxide  &  sulphide 
x  ,n 


R 

x,n 


(ki  Cel,OH,n  +  ki+l 


C  .  _  )  X 

el,S,n  n 


i.e. 


(A. 2-3) 
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where 

k_^  is  the  hydroxide  kinetic  constant  and  is  a  function 
of  temperature,  chemicals  and  wood  component  concentrations, 


^i+1  sulphide  kinetic  constant  and  is  related  in  a 

similar  manner  as  k., 

1 

Cgi  oh  n  t^ie  entrapped  liquor  hydroxide  concentration 
in  section  n, 

and  C  ,  is  the  sulphide  concentration  in  the  entrapped 

cljO jU 

liquor. 


Entrapped  liquor  component  balances  : 

The  material  balance  for  liquor  component  Y  in  the 

entrapped  liquor  is  : 

V  .  dC  ,  .  „  ^ 

_  el,n-l  el,y,n-l  el,n  el,y,n 


dt 


+  k  a  (  C-.  -C  .  ) 

m  fl,y,n  el,y,n 


(a.R.  +  a.R  +a_R  +a.R  . 

1  l,y,n  2  c,y,n  3  m,y,n  4  x,y,n) 


(A. 2-4) 

where  the  first  two  terms  are  convection  flow  terms  with  Q 

el  ,n 

being  the  entrapped  liquor  flow  rate  from  section  n,  ^ 

is  the  entrapped  liquor  concentration  for  the  component  y  in  the 
section.  The  a’s  are  the  stoichiometric  coefficients  for  the 


four  wood  components  and  liquor  chemicals  and  are  used  to  calculate 
the  chemical  usage  while  the  R's  are  the  reaction  terms  (e.g.  R 

j- » y » ti 

is  the  rate  of  reaction  of  lignin  with  chemical  y  in  section  n) . 


. 


' 


i:  . 


187 


Free  liquor  component  balances  : 

The  free  liquor  balances  for  hydroxide,  sulphide  and 
reaction  products  are  represented  by  the  general  component 
equation  : 


V - _  dC _  -  Q^-i  i^.pi  i  —  *i  C-_ 

fl,n  fl,y,n  fl,n-l  fl,y,n-l  fl,n  fl,y,n 


dt 


-  k  a  (C_  -C  .  ) 

m  fl,y,n  el,y,n 


(A. 2-5) 


these  terms  are  similar  to  those  in  the  entrapped  liquor  equation 

(A. 2-4)  with  a  convection  and  mass  transfer  term.  The  k  a  term  is 

m 

the  overall  mass  transfer  coefficient  times  area  at  the  section  temperature 


Energy  balances  : 

The  energy  balance  equations  for  wood  and  entrapped 
liquor;  and  free  liquor  follow  : 


m  o  J  rn 

t,wel,n  p  n=(w  -C  T  .  -  w  C  T)+(w-  .C  _T  n 
_  n-1  p,w  n-1  n  p,w  n  el,n-l  p,l  n-1 

dt 


-  w  ,  C  T  )  +  ZR.A/?.  +  ha  (T  -T  ) 
el ,n  p ,1  n  l  l  f ,n  n 

(A. 2-6) 


in  this  equation 

C  and  C  ,  are  specific  heats  of  wood  and  liquor 
p,w  p,l 

respectively  and  C  is  the  average  for  the  two, 

P 

w  -  is  the  entrapped  liquor  mass  flow  rate, 
el  ,n 

is  the  heat  of  reaction  for  wood  component  i, 
is  the  reaction  term, 

and  ha  is  the  heat  transfer  coefficient  times  area  between 
the  wood  and  free  liquor. 
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The  free  liquor  energy  balance  is  : 

T  7  r  JT 

fl,n  p,l  f,n  =  (w^.-  « C  _T  -  w  C  -T  ) 

**  -  f 1 ,n-l  p,l  f,n-l  f 1 ,n  p,l  f,n 

dt 


-  ha  (T,  -T  )  +  H  -| 
f,n  n  wall 


(A. 2-7) 

the  free  liquor  equation  has  similar  terms  to  the  earlier  equation 


(A. 2-6)  except  for  the  heat  conduction  from  the  walls  term, 
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A. 3  Nomenclature 

Specific 


a 

interfacial  area  between  wood  and  free  liquor 

C 

concentration  of  wood  component 

C 

0 

concentration  of  wood  component  at  the  top  of 

the  digester 

C 

n,o 

C  1 

el,y,n 

C  r: 

fl»y*n 

c  1 
pi 

c 

pw 

wood  concentration  in  section  n  at  time  zero 

concentration  of  component  y  in  entrapped  liquor 

concentration  of  component  y  in  free  liquor 

specific  heat  of  liquor 

specific  heat  of  wood  chips 

c 

p 

D 

average  specific  heat  of  wood  and  entrapped  liquor 

diffusion  coefficient 

h 

overall  heat  transfer  coefficient 

H 

total  height  of  digester 

^wall 

heat  transfer  through  walls  of  digester 

HFAC 

H-Factor,  a  measure  of  cooking  history 

k. 

l 

ki+l 

k 

m 

reaction  rate  constant  for  hydroxide 

reaction  rate  constant  for  sulphide 

overall  mass  transfer  coefficient 

mt ,wel 

Q 

total  mass  of  wood  and  entrapped  liquor  in  section 

volumetric  flow  rate 

Qel 

Qfl 

r 

volumetric  flow  rate  of  entrapped  liquor 

volumetric  flow  rate  of  free  liquor 

radius  of  section,  radial  component  in 

cylindrical  coordinates 

R,  R. 

*  i 

total  reaction  term  for  wood  component 

R 

x,n 

R . 
i 

total  reaction  term  for  wood  component  x  in  n 

heat  of  reaction 

I 


' 

. 

>• 


Specific 


w 

w 

w 


el 

fl 


X 


z 


time  variable 

wood  and  entrapped  liquor  temperature 
free  liquor  temperature 
velocity  component 
volume  of  section 

volume  of  entrapped  liquor  in  section 

volume  of  free  liquor  in  section 

wood  mass  flow  rate 

mass  flow  rate  of  entrapped  liquor 

mass  flow  rate  of  free  liquor 

concentration  of  wood  component  x 

position  vector,  axial  distance,  cylindrical 

coordinates 


Greek  Letters 


a 


1 


P 

6 


stoichiometric  coefficients  for  wood  component 
and  entrapped  liquor  chemicals 
density  of  wood  chips 

angular  component  in  cylindrical  coordinates 


Subscripts 


el 

fl 

n 

r 

e 

z 


entrapped  liquor  variable 
free  liquor  variable 
variable  from  section  n 

radial  component  in  cylindrical  coordinates 
angular  component  in  cylindrical  coordinates 
axial  component  in  cylindrical  coordinates 
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APPENDIX  B 


MODEL  TESTING  RESULTS 

This  appendix  is  related  to  the  description  of  the 
model  that  was  presented  in  Chapter  Four.  Results  from 
testing  of  integration  routines  are  presented  in  the  first 
section  and  testing  of  the  open  loop  model  and  the  kinetic 
model  in  the  other  two. 

B .  1  Testing  of  Integration  Routines 

The  four  integration  routines  tested  for  subsequent 
use  in  the  digester  model  were  variations  of  the  fourth 
order  Runge  Kutta  t^pe  and  a  predictor  corrector  method. 

The  simple  Runge  Kutta  routine  was  programmed  in  Fortran 
while  the  other  three  were  taken  from  two  subroutine 
packages  (39,  72).  The  routines  Runge  Kutta  Gill  ( EKGS )  and 
Hamming's  Predictor  Corrector  (HPCG)  were  taken  from  the  IBM 
Scientific  Subroutine  Package  (SSP)  and  the  Runge  Kutta 
Merson  (DASCRU)  from  the  IMSL  package.  The  routines  were 
used  to  integrate  the  equations  ten  times  sc  that  a 
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reasonable  values  for  CPU  time  could  be  obtained. 

The  first  set  of  equations  used  in  testing  the 

routines  were  linear  and  the  analytical  solution  of  which 

was  in  exponentials  . 

The  equation  set  was  : 

dY  (l)  =  _  1J 

dt  “  Y(2) 

<11  (2)  = 

dt  Y(1)  (B.1.1) 

and  the  conditions  at  time  zero  were  Y(1)  =  Y (2)  =  1.0.  The 
analytical  solution  for  this  problem  is  : 

Y  ( 1 )  =  exp(-t) 

and  Y  (2)  =  exp  (  +  t)  { B -  1-2) 

Table  B. 1  shows  the  computing  time  used  in  CPU 


seconds  on 

an  Amdahl 

47C/V6  computer  and  the 

solutions 

obtained  from  the  integration  routines. 

The  second 

set  used  was  a  slightly 

stiff  set  - 

eq  uations , 

and  they 

are  as  follows: 

dY 

dt 

(1) 

-0. 1*Y (1)  -  49 . 9*Y  (2) 

dY 

dt 

(2) 

-50.  0*  Y  (2) 

dY 

(3) 

70. 0*  Y  (2)  -  1 2  0. 0* Y  (3) 

(3. 1.2) 

dt 

. 
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The  initial  conditions  in  this  case  were,  Y  (1)  =  Y (3)  =  2.0 
and  Y (2)  =  1.0.  This  equation  can  be  solved  analytically 
also  and  the  integration  results  for  this  test  are  given  in 
Table  B.2.  Based  on  these  two  tests  the  Runge  Kutta  Gill 
and  Predictor  Corrector  routines  were  definitely  inferior  to 
the  simple  Runge  Kutta  and  Runge  Kutta  Merson  methods  based 
on  accuracy  and  computing  time  usage. 

The  third  test  was  on  the  digester  model  equations 
where  the  final  solution  was  known.  This  was  one  of  the 
trivial  tests  with  no  reaction  and  also  no  heat  and  mass 
transfer  between  the  wood  and  liquor.  The  Runge  Kutta  Gill 
and  Predictor  Corrector  routines  were  not  tested  for  this 
case.  Two  profiles  were  calculated  and  results  from  both 
methods  were  identical.  The  CPD  time  for  the  Runge  Kutta 
was  1.750  seconds  and  that  for  the  Merson  method  3.235 
seconds.  This  showed  that  the  simple  Runge  Kutta  routine 
was  performing  much  better  than  the  other  routines  when  the 
problem  was  simple  and  the  nonlinearities  not  very  severe. 
The  actual  digester  equations  are  extremely  nonlinear 
especially  in  the  vicinity  of  the  heater  sections  where  the 
pulp  temperature  is  raised  by  at  least  thirty  degrees 
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Celsius-  The  Runge  Kutta  routine  was  not  able  to  outperform 
the  Merson  routine  in  this  case  and  very  small  step  sizes 
had  to  be  used  to  obtain  reasonable  results-  The  Merson 
routine  was  for  this  very  reason  used  to  solve  the  nonlinear 
equations  in  the  batch  and  continuous  digester  models. 
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0. 1 92870E-2 1  0, 3720  10E-43 
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• 

B.  2  OjJen  Loop  Model  Testing 

The  program  subroutines  were  tested  in  parts  and 
the  test  and  results  were  as  follows  : 

Table  E.3  Open  Loop  Program  Testing 


PROGRAM  SEGMENT 

TEST 

RESULT 

Digester 

Configuration 

Calculation 

Volumes  and  Areas 

calculated 

manually 

Results  agreed 
perfectly 

Evaluate  liquor 
density  A  its 
variation  with 
tempe  ratur e 

Density  values 
calculated  at 
different 
temperatures 

Values  were  within 
0- 1  %  of  those  in 
the  International 
Critical  Tables 

Sodium  Sulphide 
hydrolysis  routine 

Varying  NaOH  A 
Concentrations 
Temperature 

Na2  S 
and 

Results  checked 
with  Martin's  (o4) 
data.  Errors  were 

large  (Max  12%)  at 
lew  temperatures  A 
high  sulphidities 


con  tinued 
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Table  E.3  continued 


PROGRAM  SEGMENT 

TEST 

RESULT 

Sodium  Sulphide 
available  to 
degrade 
carbohydrates 

Print  out  of 
values  at  different 
process  conditions 

Subroutine  worked 
well.  Part  of  the 
batch  model 

Residual  reaction 
phase  breakpoint 

Check  subroutine 
working 

ii  it 

Reaction  velocity 
constants 

Kinetic  constants 
at  various  process 
conditions 

Values  agreed 
fairly  well  with 
Stillman  (81)  data 

Mass  transfer 

Coef f icient 
calculations 

Checks  at  various 
temperatures 

Data  from  Smith  & 
Williams  (76)  . 

Perfect  mixing 
calculations  for 
external  flows 

Hand  calculations 
for  material  and 
energy  balances 

Calculations 

match 

Calculation  of 
Compaction  values 
as  a  function  of 
position  &  history 

Values  for  position 
down  the  digester 
&  with  various 
cooking  histories 

Agree  well  with 
Hamilton^  (29) 
data 

Digester  entering 
conditions  at  the 
top 

Check  working  with 
hand  calculations 

Results  match 

continued 


. 


« 

* 

■ 


199 


Table  E.3  continued 


PROGRAM  SEGMENT 

TEST 

RESULT 

Differential 

Test 

Equation  subroutine 

on 

sample 

Results 

Runge  Kutta  Merson 

equations 

are 

method  for  solution 

and 

as 

of  equation  set 

for 

expected. 

trivial 

i.e. 

Routines  for  solving 

tests. 

const  ant 

cocurrent  &  counter- 

wood 

current  flow  problem 

i.e.  no  reaction  & 

component 

mass  transfer,  no 

profiles  & 

Mainline,  driver 

heat  transfer 

entrapped 

program 

liq  uor 

concentrations 

< 
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B.3  Kinetic  Model  Testing  -  Batch  Model  Working 

The  Kinetic  Model  that  was  described  in  Chapter 
Three  was  tested  by  using  a  batch  digester  nodel  that  was 
written  specifically  for  that  purpose.  The  batch  model  had 
the  same  type  of  thirteen  nonlinear  differential  equations 
as  in  the  continuous  model  except  that  the  flow  terms  were 
deleted.  The  results  from  this  model  were  tested  against 
cooking  data  collected  by  Daleski  (19)  and  Smith  and 
Williams  (88)  and  the  kinetic  parameters  were  tuned 
accordingly  to  give  the  best  possible  fit.  The  computer 
model  is  such  that  if  better  kinetic  data  for  a  particular 
species  were  available  it  could  be  read  in  and  used  for 
simulations.  The  constants  given  in  Table  2  are  average 
values  for  hardwood  and  softwood  species  and  were  used  to 
generate  the  plots  that  are  shown  in  Figures  B.  1  to  B-5. 


. 
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Figure  B.l 


Batch  Cook,  Yield  at  195  °C  (Daleski  (19)) 


.  LIGNIN.  V.  WOOD 

3..7S  7.50  1.1.25  15.00  18.75  22.50  26.25  30.00 
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Batch  Cook,  Lignin  at  195  °C  (Daleski  (19)) 


Figure  B.2  : 


.  CELLULOSE.  WOOD 

5.63  11.25  16.88  22.50  28.13  33.75  39.38  45.00 


203 


o 

o 


4U0 


6*.  25 


l1. 00 


^.00  0.75  1  .50  2.25 

COOKING  TIME 


3.00  3.75 

MINUTES 


Figure  B.3 


Batch  Cook,  Cellulose  at  195  °c 


(Daleski  (19)) 


,  MRNNflN »  Z.  WOOD 

1.88  3.75  5.63  7.50  9.38  11.25  13.13  15.00 
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*^.00 


l‘.BO  2*25  3*.  00  3V75  4 '.50  s‘-25  6.00 

COOKING-  TINE,  NINUTES 


Eatch  Cook,  Mannan  at.  195  °C  (Daleski  (19)) 


Figure  B.4  : 


XYLAN .  V.  WOOD 

2.00  3.00  4..00  5.00  -  6..00  7.00  8.00 
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COOKINO  TIME.  MINUTES 


Figure  B.5  :  Batch  Cook,  Xylan  at  195  °C  (Daleski  (19)) 


APPENDIX  C 


WORKING  OF  OPEN  LOOP  PROGRAM 


This  appendix  contains  the  flow  diagram  of  the  main 
computer  pregram  and  a  set  cf  sample  data  with  the  standard 
program  output-  The  dynamic  and  steady  state  digester 
models  are  a  part  of  one  program  which  is  atout  three 
thousand  lines  long  and  the  FORTRAN  G  compiled  version  on 
the  Amdahl  470/V6  computer  occupied  a  256K  bytes-  The  run 
time  to  generate  the  profiles  in  this  appendix  was  2-5  CPU 
seconds  when  the  program  was  already  loaded  in  memory  and 
3.0  CPU  seconds  when  it  was  not.  The  batch  program  is  only 
1500  lines  long  and  its  compiled  version  occupied  atout  100K 
bytes. 


C-  1 


expanded 
shown  in 
complete 
manual . 


Open  Loop  Model  Flow  Diagram 

The  flow  diagram  in  this  appendix.  Figure  C. 1  is  an 
version  of  the  simplified  flow  diagram  that  was 
Figure  4.1.  The  more  detailed  flow  diagram  with 
program  documentation  is  given  in  the  programmers 
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Figure  C.l 


Open  Loop  Model  Flow  Diagram 
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Figure  C.l  ..  continued 
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C . 2  Ogen  Loop  Sample  Data  and  Printout 


The  data  used  to  produce  the  results  for  the  sample 


computer  run  in  the  this  section  is  presented  here. 


Setup  of  Program  Flags 


Wood  Type 
Wood  Species 
Units 

Initial  Profile  Type 
History  Required  ? 

Heat  Balance  ? 

Secondary  Coox? 

Equation  Testing? 

Changes  in  Operating 
Conditions 
Soda  Cook? 

Rate  Constants  Readin 
Reaction? 

Washing  Calculations? 

Profile  Printing 
Frequency 

Profile  Save  Frequency 
Plotting? 

Debugging  Messages? 

Stopping  Criterion? 

**  DATA  ** 

Blow  Flow 

Pulp  Consistency 

Chip  flow 

Chipmeter  Capacity 
Chip  Moisture 
Chip  Critical  Moisture 
Wood  Specific  Area 


Hardwood 

Beech 

System  Internationale  (S.I.) 

Startup  from  scratch 

No 

No 

No 

No 


No 

No 

Yes 

Yes 

Yes 


All 

Every 

Fourth 

Yes 

No 

After 

TOT  CEO  Profiles 

54. 

0 

cu 

m/min 

12. 

5 

% 

63500 

.  0 

kg/hr 

46 

kg/ re 

volution 

50. 

0 

% 

60. 

0 

% 

155 

.  0 

sq 

m/cu  m 

' 


. 

'  -- 
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Dry  Wood  Density 
Wood  Substance  Density 
Lignin  Content 
Cellulose  Content 
Mannan  Content 
Xylan  Content 
Acetyls 
Uronic  Acid 
Extractives 

Percent  Extractive  Lignin 
Percent  Extractive  Mannan 
Percent  Extractive  Xylan 
White  Liquor  Flow 
White  Liquor  Alkali 
White  Liquor  Sulphidity 
White  Liquor  Solids 
White  Liquor  Temperature 
Black  Liquor  Flow 
Black  Liquor  Alkali 
Black  Liquor  Sulphidity 
Black  Liquor  Solids 
Black  Liquor  Temperature 
Upper  Beater  Flow 
Lower  Heater  Flow 
Quench  Flow 
Wash  Heater  Flow 
Upper  Heater  Temperature 
Lower  Heater  Temperature 
Quench  Temperature 
Wash  Heater  Temperature 
Extraction  Flow 


184.2  kg/cu  m 
1653.0  kg/cu  m 
26.6  kg/IOOkg  ODW 
39.0  kg/IOOkg  ODW 

16.2  kg/IOCkg  ODW 

9.5  kg/IOOkg  ODW 
0.95  kg/IOOkg  ODW 
3.85  kg/IOOkg  ODW 
3.90  kg/IOOkg  ODW 
15.0  % 

55.0  % 

15.0  % 

12.75  cu  m/min 

140.5  kg/cu  m  as  Na20 
33.0  % 

0.0  % 

87.0  C 

0.0  cu  m/min 
8. 0  kg/cu  m 
15.0  % 

3.9  % 

85.0  C 

3  0.5  cu  m/rain 
28.0  cu  m/min 

Calculated  Latter  in  Program 

9.0  cu  m/min 

150.0  C 

155.0  C 

98.0  C 

94.0  C 

35.0  cu  m/min 


Quench  Flow/Extraction 
Ratio 

Cold  Blow  Alkali 
Cold  Blow  Sulphidity 
Cold  Elow  Solids 
Cold  Elow  Temperature 
Kinetic  Parameters 


Flow 

0.40 

0.0  kg/cu  m 

0.0  % 

0.0  % 

45-0  C 

Same  as  in  Table  2. 


Heat  Transfer  Coefficients 
Internal  Pipe  and  Liquor  0.235  W/sq  m  -  K 
Digester  Walls  0.0002  W/sq  m  -  K 

Free  Liquor  and  Wood  4.72  W/sq  m  -  K 


. 


T 
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Heat 


Heat 

Heat 


of  Reaction 

580.0  J/kg 

of  wood 

fic  Heat  of  Wood 

2.878  J/kg 

-  K 

fic  Heat  of  Liquor 

5.  188  J/kg 

-  K 

fic  Heat  of  Ice 

2.090  J/kg 

-  K 

of  Fusion  of  Ice 

93.0  J/kg 

of  Condensation  of 

Steam  2188.0  . 

J/kg 

rature  of  Chips 

30.0  C 

rature  of  Steaming 

Vessel  125.0  C 

nt  Temperature 

29.0  C 

reacted 


Upper 

Lower 


Cold 


ter  Dimensions 

Printed 

out  in  Sample  Outpu 

Number  of  Sect 

ions  per 

Zone 

gnation 

3 

Heater 

2 

Heater 

Always 

1 

ng 

6 

h 

Always 

1 

ction 

Always 

1 

ng 

7+1  for 

Wash  heater  screen 

Elow 

2 

ical  Step 

Size 

5.0  min 

les  to  be 

Calculated 

2 

Checking 

constants  fo 

r  Variables 

Profile 

Plotted  ? 

No 

Every 

Number  of  Variables  Plotted  6 
Total  Number  of  Plots  1 
Profile  Number  of  Plot  1 


Variables  to  be  Plotted 
Lignin 
Cellulose 
Mannan 
Xylan 

Entrapped  Hydroxide 
Chip  Temperature 


' 


■ 
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C .3  Documentation  Produced  by  the  Computer  Program 

The  sample  computer  run  follows  and  consists  of  tne 
digester  configuration,  volumes  and  areas  in  each  section, 
initial  digester  entering  flows,  the  initial  digester 
profile  and  the  dynamic  profiles.  The  plots  and  other 
output  from  a  similar  run  were  shown  in  chapter  five. 

Details  are  given  in  the  programmer’s  manual. 


^  • 


Table  C.l  :  Open  Loop  Sample  Data  and  Run 
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«  •  • 


continued 


TOP  OF 

cot;  e 
4.  10 
4.60 


DP(1) 

37.50 

ZONE 

IKPG 

OPHT 


LBHT 


CKNG 


QUEN 


EXTB 


WASH 


♦  »*»  CONFIGURATION  OF  THE  DIGESTER  **M 


TOTAL  HEIGHT  OP  DIGESTER,  3 
70.  30 

TOTAL  VOLUME  CU.  H 
1265.86 


ZONE  HEIGHTS,  M 


IH PHEG  NATION  14. eO 
OPPEP  HEATING  3.25 
LOWER  HEATING  1.90 
COOKING  18.50 
QUENCH  2.35 
EXTRACTION  2-60 
WASHING  23.00 
COLD  BLOB  3. SO 


INSIDE  DIAMETER  AND  HEIGHT  OF  DIGESTER  ZONES,  H 


BOTTOM 

D{3) 

DJ4) 

D (5)  D (6 ) 

D{7) 

D(8) 

D(9) 

0?  CONE 

4.85 

4.  05 

4.95 

5.00  5.10 

4.  85 

4.80 

4.75 

10.00 

3.25 

1.90 

18.50  2.35 

2.60 

23-0  0 

3-90 

OUTSIDE  DIAMETERS  OF  PIPE 

BUNDLE, 

CM 

DP  (2) 

DP(3) 

DP  (4) 

DP  (5) 

DP(6) 

DP  (7) 

DP  (8) 

55.75 

50.00 

5  0.00 

45.00 

4  0.00 

40.00 

40.00 

*♦*  SECTION  CONFIGURATION  *** 


SECTION  AN? 

VOLUME 

PIPE  AREA 

WALL  AREA 

SCREEN  NO. 

CU.  H 

SQ.  H 

SQ.  B 

1 

32-61 

2.64 

29.56 

2 

25.03 

4.10 

20.91 

3 

181.68 

23.29 

152. 37 

4 

29-98 

4.53 

25.  14 

SCREEN  1 

5 

29-07 

4.40 

24.38 

6 

31.25 

4.53 

25.66 

SCREEN  2 

7 

59.71 

7.99 

48.43 

8 

59.71 

7-99 

48.43 

9 

59.71 

7.99 

48.43 

10 

59.71 

7.99 

48.  43 

11 

59-71 

7.99 

48.43 

12 

59.71 

7.99 

48.  43 

13 

47.71 

2.95 

37.65 

SCREEN  3 

14 

47.71 

3.27 

39.62 

SCREEN  4 

15 

56.09 

3-92 

47.07 

16 

56.09 

3.92 

47.07 

17 

56.09 

3.92 

47.07 

18 

56.09 

3.92 

47.07 

19 

56.09 

3.92 

47.07 

20 

56.09 

3.92 

47.0  7 

21 

56.  09 

3.92 

47.07 

22 

SCREEN  5 

20.81 

0.0 

17.34 

23 

34.56 

0.0 

29.  10 

24 

34.56 

0.0  , 

29.10 

CBLO 
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